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L  INTRODUCTION 

This  thesis  has  two  primary  goals.  The  first  goal  is  to  derive  and  document  the 
solution  to  the  three-dimensional  inhomogeneous  Helmholtz  equation  for  a  ffee-space 
propagation  problem  when  the  square  of  the  index  of  refraction  is  a  linear  function  of 
depth.  Given  this  type  of  ocean  medium,  the  inhomogeneous  Helmholtz  equation  has  an 
exact  solution  in  terms  of  Airy  functions.  This  exact  "Airy  function"  solution  will  be 
incorporated  into  the  computer  program  Linear  Space- Variant  Ocean  (LSVOCN)  as  an 
additional  ocean-medium  transfer  function.  The  second  goal  of  this  thesis  is  to  further 
test  the  Recursive  Ray  Acoustics  (RRA)  Algorithm  by  comparing  the  magnitude  and 
phase  of  the  acoustic  pressure  calculated  along  a  ray  path  by  the  RRA  Algorithm  with  the 
magnitude  and  phase  calculated  by  LSVOCN. 

The  first  step  in  this  process  requires  a  careful  derivation  of  the  solution  to  the 
three-dimensional  inhomogeneous  Helmholtz  equation  for  a  free-space  propagation 
problem  when  the  square  of  the  index  of  refraction  is  a  linear  function  of  depth.  This 
derivation  will  be  discussed  in  detail  in  Chapter  II  which  is  divided  into  three  main 
sections.  Section  A  defines  the  three-dimensional  inhomogeneous  Helmholtz  equation  to 
be  solved.  (By  using  a  zeroth-order  Hankel  transform,  the  solution  process  is  simplified 
to  solving  a  one-dimensional  inhomogeneous  ordinary  differential  equation  vice  a 
three-dimensional  inhomogeneous  Helmholtz  equation).  Section  B  defines  the 
mathematical  formula  used  to  model  the  square  of  the  index  of  refraction  as  a  linear 
function  of  depth.  Section  C  documents  the  three  steps  that  lead  to  the  final  solution  of 
the  inhomogeneous  Helmholtz  equation. 

First,  the  description  of  the  sound-speed  profile  is  defined.  This  description 
allows  us  to  solve  for  the  two  constants  in  the  equation  for  the  square  of  the  index  of 
refraction  as  a  linear  function  of  depth.  Second,  the  homogeneous  depth  equation  is 
manipulated  into  the  form  of  Airy’s  differential  equation  and  then  solved  in  terms  of  Airy 
and  Bairy  functions.  Third,  the  solution  of  the  inhomogeneous  depth  equation  is  obtained 
by  applying  the  appropriate  boundary  conditions  at  the  source. 
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The  solution  to  the  inhomogeneous  Helmholtz  equation  is  divided  into  two  cases. 
In  case  one,  the  square  of  the  index  of  refraction  has  a  negative  gradient.  This  case  is 
further  subdivided  into  a  solution  for  the  receiver  depth  equal  to  or  above  the  source 
depth,  and  a  solution  for  the  receiver  depth  below  the  source  depth.  In  case  two,  the 
square  of  the  index  of  refraction  has  a  positive  gradient,  and  again  a  solution  for  the 
receiver  depth  equal  to  or  above  the  source  depth,  and  a  solution  for  the  receiver  depth 
below  the  source  depth  are  found. 

Chapter  III  is  divided  into  three  sections.  Section  A  discusses  the  numerical 
techniques  used  in  the  process  of  converting  the  Airy  function  solution  into  working 
FORTRAN  computer  code  used  by  LSVOCN,  and  defines  the  Green's  function  used  to 
test  LSVOCN.  Before  testing  of  the  RRA  Algorithm  began,  a  comparison  was  made  of 
the  output  of  LSVOCN,  using  the  above  derived  Airy  function  solution,  with  a  test  case, 
in  an  effort  to  validate  the  FORTRAN  computer  code.  The  simplest  test  case  considered 
compared  the  magnitude  and  phase  of  the  time-independent  free-space  Green's  function 
for  an  isospeed  medium  with  the  output  from  LSVOCN. 

Section  B  of  Chapter  III  discusses  the  three  different  Green’s  function  test  cases 
that  were  used  to  test  LSVOCN.  These  are:  (1)  Green's  function  test  case  with  the  source 
sound-speed  varied  such  that  the  modeled  ocean  medium  approached  an  isospeed 
medium  in  order  to  compare  the  magnitude  and  phase  of  the  Green's  function  with 
LSVOCN,  (2)  Green's  function  test  case  with  the  horizontal  range  varied  while 
maintaining  the  receiver  at  the  source  depth,  thus  simulating  an  isospeed  ocean  medium 
for  comparison  of  the  Green's  function  magnitude  and  phase  with  LSVOCN,  and  (3) 
Green's  function  test  case  with  the  receiver  depth  varied  above  and  below  the  source 
depth  by  one-half  meter  allowing  for  an  approximate  comparison  of  the  Green's  function 
magnitude  and  phase  with  LSVOCN.  In  each  test  case,  the  change  in  sound  speed 
between  the  source  and  receiver  was  small  enough  to  approximate  an  isospeed  ocean 
medium. 

Section  C  of  Chapter  III  focuses  on  the  validation  of  the  predicted  values  for  the 
magnitude  and  phase  of  the  acoustic  pressure  along  a  ray  path  given  by  the  RRA 
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Algorithm.  This  is  done  by  propagating  sound  rays  at  various  launch  angles  using 
negative  and  positive  gradients  for  the  index  of  refraction  profile.  The  horizontal  range 
and  depth  coordinates  of  selected  points  along  these  ray  paths  have  been  recorded  and 
will  become  the  receiver  horizontal  range  and  depth  coordinates  used  later  as  inputs  to 
LSVOCN.  The  magnitude  and  phase  values  calculated  by  LSVOCN  are  then  compared 
to  the  corresponding  magnitude  and  phase  values  obtained  from  the  RRA  Algorithm  for 
validation  purposes. 

Chapter  IV  summarizes  the  numerical  problems  encountered  during  the  validation 
process,  states  conclusions  as  well  as  recommendations  for  future  research. 
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n.  AIRY  FUNCTION  SOLUTION  OF  THE  INHOMOGENEOUS 
HELMHOLTZ  EQUATION 

A.  THE  INHOMOGENEOUS  HELMHOLTZ  WAVE  EQUATION 

In  this  chapter  we  shall  derive  a  solution  for  the  following  inhomogeneous 

Helmholtz  equation  in  cylindrical  coordinates  (r,  <j>,y),  shown  in  Fig.  1: 

dr2  ‘P/fojO +  r  lrVf  (w) +  ^2  <P/(r^) + *2(y)<pAr,y)  = 

(2.1) 

where  the  velocity  potential  yf(r,y)  is  axisymmetric  (i.e.,  independent  of  the  azimuthal 

angle  <p ), 

(2.2) 

is  the  depth-dependent  wave  number,  c(y)  is  the  depth-dependent  speed  of  sound,  and  the 
impulse  functions  on  the  right-hand  side  of  Eq.(2.1)  represent  a  unit  amplitude, 
omnidirectional  point  source  at  r  -  0  and  7  =y0. 

Figure  1.  Illustration  of  problem  geometry. 
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The  solution  tyf(r,y)  of  Eq.(2.1)  can  be  expressed  as  an  inverse  zeroth-order 
Hankel  transform  as  follows,  see  [Ref.  2,  Appendix  3C]: 

c Mr,y)=Hj{<bf{kr,y)}  =  \~Q>j{kr,y)J0(krr)krdkr  (2.3) 

where 

,y)  =H0{<p f{r,y)}  =  Jq  <p/(r,y)Jo(krr)rdr  (2.4) 


is  the  forward  zeroth-order  Hankel  transform  (also  known  as  the  Fourier-Bessel 
transform)  of  (fy(r,y).  Therefore,  if  we  can  find  <&f{kr,y),  then  substituting  <bf(kr,y)  into 
Eq.(2.3)  will  yield  a  solution  of  Eq.(2.1). 


and 


Since 

+  T^M,y))  =  -k2r<&f(Ky) 


(2.5) 

(2.6) 


taking  the  zeroth-order  Hankel  transform  of  Eq.(2.1)  yields 


^®Akr,y)  +  k2y(y)<bf{kr,y)  =  5^2^o) 


(2.7) 


where 


k2yiy)  =  k2(y)-k2.  (2  j 

In  Eq.(2.8),  ky{y)  is  the  depth-dependent  propagation-vector  component  in  the  Y 
direction  and  kris  the  constant  propagation -vector  component  in  the  horizontal  range  r 
direction  (see  Fig.  2).  Therefore,  now  we  only  have  to  solve  the  ordinary  differential 
equation  given  by  Eq.(2.7)  instead  of  the  partial  differential  equation  given  by  Eq.(2. 1 ). 
Substituting  the  solution  <t>f{kr,y)  into  Eq.(2.3)  yields  the  solution  yf{r.y)  to  Eq.(2.1). 
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Figure  2.  Illustration  of  the  propagation-vector  components  kr  and  ky  iy) 
and  their  relationship  to  the  wave  number  k(y). 


B.  PROBLEM  CONSTRAINTS 

The  solution  of  Eq.(2.1)  that  shall  be  derived  is  only  applicable  to  ffee-space 
propagation  when  the  square  of  the  index  of  refraction  is  a  linear  function  of  depth  y.  The 
index  of  refraction  is  defined  as  follows: 

<2-9> 

where  c0  =  c(y0)  is  the  speed  of  sound  (m/s)  at  the  source  depth  j>0  and  c(y )  is  the  speed 
of  sound  (m/s)  at  depth  y.  If  we  let  k0  =  k(y0),  then  evaluating  Eq.(2.2)  aty  =  y0  yields 


(2.10) 


Therefore,  we  can  write  that 


*0)  =  k0n(y)  = 


Inf 

c(y) 


(2.11) 


As  a  result,  Eq.(2.8)  can  be  rewritten  as 
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(2.12) 


by  (y)  =  k\n7  (y)-k*. 

To  ensure  that  the  square  of  the  index  of  refraction  is  a  linear  function  of  depth,  we 
define  the  following  equation: 

n2(y)=aiy  +  a0.  (2.13) 

C.  PROBLEM  SOLUTION 

Using  the  results  from  the  last  two  sections,  specifically  Eq.  (2.7)  and  (2.12),  we 
are  now  ready  to  solve  the  inhomogeneous  Helmholtz  equation.  The  solution  can  be 
broken  into  two  parts.  First,  determine  the  solution  of  Eq.(2.7)  aty  =y0 ,  i.e.,  solve  the 
homogeneous  equation.  Solving  the  homogeneous  equation  will  result  in  four  unknowns. 
In  order  to  determine  these  four  unknowns,  which  is  the  second  part  of  the  problem,  we 
apply  the  boundary  conditions  above  and  below  the  source  and  the  boundary  conditions 
at  the  source,  i.e.,  aty  =  y0.  As  shown  in  [Ref.  2,  Appendix  3C],  the  two  boundary 
conditions  that  must  be  satisfied  at  the  source  are  continuity  of  acoustic  pressure  and 
discontinuity  in  they  component  of  the  acoustic  fluid  velocity  vector.  With  this  brief 
introduction  in  mind,  let  us  now  proceed  with  the  solution  of  the  inhomogeneous 
Helmholtz  equation. 

L  Description  of  the  Sound-Speed  Profile 

Consider  a  sound-speed  profile  where  the  square  of  the  index  of  refraction  is  a 
linear  function  of  depth,  [see  Eq.(2.13)]: 

n1(y)=a\y  +  ao.  (2.14) 

With  the  above  equation  in  mind,  we  must  now  solve  for  the  two  unknown  constants, 
a0  and  av  in  order  to  solve  Eq(2.7).  Substituting  Eq.(2.9)  into  Eq.(2.14)  yields 

cl 

7^)-aiy  +  a°-  (2-15) 
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Rewriting  Eq.(2. 15)  yields 


and  if  we  let 


_J_  =  £i.v  +  £o 

c2(y)  cf  cV 


/  _a  i 

ai=~I 


(2.16) 


(2.17) 


/  _ 

a0  ”  2  » 

Co 


(2.18) 


i=a^  +  a°- 


c2(y) 


(2.19) 


Solving  for  the  speed  of  sound  yields 


c(y)  = 


a\y+a'0  ^/uriy+ao 


(2.20) 


Equation  (2.20)  will  allow  the  sound-speed  profile  to  be  calculated,  once  the  unknown 
constants  have  been  determined. 

Evaluating  Eq.(2.19)  aty  =  0  and  at y  =y0  ,  i.e.,  the  source  depth,  yields 


/  _  1 
°  c2(0) 


(2.21) 


r  r  —  aiy° + c'o  > 


c2Oo) 


(2.22) 


respectively.  Substituting  Eq.(2.21)  into  Eq.(2.22)  and  solving  for  a\  yields 


/  1  [  J _ L 

1  c2(( 


(2.23) 
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Thus,  using  Eqs.(2.21),  (2.23),  (2.17),  and  (2.18),  a0  and  a ,  can  be  determined. 

For  our  problem,  we  will  use  a  source  depth  of  v0  =  150  meters.  In  order  to 
define  a  sound-speed  profile,  we  specify  the  sound  speed  (m/s)  at}'  =  0  to  be  c(0)  =  1500 
m/s.  For  a  positive  sound-speed  gradient  we  specify  the  sound  speed  at  the  source  depth 
to  be  c0  =  1515  m/s.  Similarly,  for  a  negative  sound-speed  gradient,  we  specify  the 
sound  speed  at  the  source  depth  to  be  c0=  1485  m/s.  These  sound-speed  profiles,  along 
with  the  corresponding  index  of  refraction  plots,  are  shown  in  Fig.  3  and  Fig.  4, 
respectively. 


Fig.  3  Positive  sound-speed  gradient. 


Fig.  4  Negative  sound-speed  gradient. 
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2.  Solution  to  the  Homogeneous  Depth  Equation 

Having  carefully  defined  the  sound-speed  profile,  we  are  now  ready  to  solve  the 
inhomogeneous  depth  dependent  Helmholtz  equation  as  given  by  Eq.  (2.7).  As  stated 
previously,  the  first  part  of  the  solution  process  is  to  solve  the  homogeneous  equation 
where  the  square  of  the  index  of  refraction  is  a  linear  function  of  depth.  The 
homogeneous  equation  is 

fys>Akr,y)+%()’)<s>A*r.y)= 0-  (2.24) 

Substituting  Eq.  (2.13)  into  Eq.  (2.12)  yields 


kj(y)  =  k%(aiy+ao)-k*, 

(2.25) 

k2(y )  =  klaxy  +  kla0  -  k?r . 

(2.26) 

Substituting  Eq.  (2.26)  into  Eq.  (2.24)  and  simplifying 

-^4>/(kr,y)  +  {axkly  +  aok20  -  kj )<$/(kr,y)  =  0 .  (2.27) 


If  we  let  oto  and  di  be  defined  as 


1 

<3 

II 

o 

a 

(2.28) 

and 

- 

oil  =a\kl , 

(2.29) 

then 

d1 

—®j(kr,y)  +  (aiy  +  a  0)®/(kr,y)  =  0. 

(2.30) 

Next,  we  must  manipulate  Eq.  (2.30)  into  the  form  of  Airy's  differential  equation. 
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(2.31) 


The  first  step  in  this  process  is  to  define  the  new  term 

_2 

C(y)  =±(«i)  3(aij  +  a0). 


Applying  the  chain  rule  yields 


and  since 


dy 


(2.32) 


(2.33) 


upon  substituting  Eq.(2.33)  into  Eq.  (2.32),  we  obtain 


(2.34) 


Applying  the  chain  rule  again  yields 


=  H^k"y) = ![*<“' : >*  (235) 


(2.36) 


and  upon  substituting  Eq.  (2.34)  into  Eq.  (2.36),  we  obtain 


&»/<*'■*>  =  +<«i)3 


^20) 


(2.37) 


Combing  Eqs.  (2.30),  (2.31),  and  (2.37)  yields 


+«*  i )  ■ 3  +  (±)(a  i ) 3 1 \{y)<S>f{kr,y )  =  0 , 


(2.38) 
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or 


jf-ftoMr.y)  ±  =  0-  (2-39) 

Upon  choosing  the  minus  sign,  we  have  Airy's  differential  equation: 

^y4>/(*„>-)-SO)4>/(*„r)  =  0.  (2.40) 

The  general  solution  of  Airy's  differential  equation  is  given  by 

Qfikr^^AAi&m+BBi&iy)]  JaryZyo  (2.41) 

*A*r>y)  =  CAi&(y)]  +  DBA C(y)]  ,  for  y  >y0  (2.42) 


where  A i  and  Bi  are  the  Airy  and  Bairy  functions,  and  A,  B,  C,  and  Dare  unknown 
constants.  Equation  (2.41)  and  Eq.(2.42)  form  the  solution  to  the  depth  dependent 
homogeneous  Helmholtz  equation.  We  have  an  equation  with  four  unknown  constants  as 
a  result  of  solving  a  second-order  differential  equation.  To  find  these  four  unknown 
constants,  we  must  apply  the  boundary  conditions  above  and  below  the  source  and  the 
boundary  conditions  at  the  source,  as  will  be  done  in  the  next  section. 

3.  Solution  to  the  Inhomogeneous  Depth  Equation 

In  this  section  we  will  determine  the  general  solution  to  the  inhomogeneous  depth 
dependent  Helmholtz  equation.  As  stated  previously,  in  order  to  solve  for  the  four 
unknown  constants  in  Eq.(2.41)  and  Eq.(2.42),  we  must  satisfy  the  boundary  conditions 
above  and  below  the  source  and  the  boundary  conditions  at  the  source.  However,  before 
we  proceed  with  this  solution,  let  us  examine  in  detail  the  behavior  of  the  functions 
yfz[£(y)]  and  /?/[C(y)]  to  gain  insight  into  the  solution  of  the  problem. 

By  choosing  the  minus  sign  in  Eq.(2.31),  we  can  write  that 

CCv)  =  — (oti  )~^(otijv  -i-oto)  (2.43) 


13 


Because 

Pyiy)  =  a.\y  +  a.0, 

Eq.(2.43)  can  be  rewritten  as 

C(y)  =  -(ai)~3*J(y). 

Furthermore,  by  referring  to  Eq.(2.12),  we  can  write  that 


(2.44) 

(2.45) 


ky(y)  =  k0[n2(y)  -  sin2(p0)] 


(2.46) 


since  kr  k0sin((3o ),  where  (i0  is  the  launch  angle,  measured  with  respect  to  the  positive  Y 
axis.  Values  of  ft  range  from  0  to  180  degrees,  [Ref.  2,  Chapter  5],  Tables  18  through 
23,  located  in  the  Appendix,  contain  values  for  n2(y),  ky2(y),  and  C,(y)  as  a  function  of 
frequency,  ft,  and  depth  y,  for  the  analysis  of  the  functions  ^/[^(y)]  and  2K[£G0].  The 

data  in  these  tables  is  further  displayed  in  graphical  form  in  Figures  12  through  1 7,  also 
located  in  the  Appendix. 

To  analyze  the  behavior  of  Ai[^(y)]  and  i?/[£(y)]  >  let  us  draw  our  attention  to  Fig. 
5  and  Fig.  6,  as  shown  below: 


Y 


Y 


The  free-space  solution  can  be  broken  into  two  cases.  Case  one:  n2  (y)  has  a  negative 
gradient,  (see  Fig.  5,  Fig.  3,  Fig.  15  through  Fig.  17,  and  Tables  21  through  23).  Case 
two:  n2(y)  has  a  positive  gradient,  (see  Fig.  6,  Fig.  4,  Fig.  12  through  Fig.  14,  and  Tables 
18  through  20).  Also,  as  can  be  seen  from  Fig.  5  and  Fig.  6  ,  each  case  must  be  analyzed 
in  two  regions,  one  above  the  source  depth  y0 ,  and  one  below  the  source  depth  y0. 

First,  let  us  analyze  the  behavior  of  the  functions  .<4i[C(y)]  and  2?/[C(y)]  for, 

Region  lib,  Case  one.  Notice  from  Fig.  5,  that  the  value  of  n2  (y)  is  less  than  one  in  this 
region.  Thus,  as  can  be  seen  from  Eq.(2.46),  ky2(y)  will  take  on  negative  values.  This  will 
result  in  values  of  £(y)  taking  on  positive  values  [see  Eq.  (2.45)].  Therefore,  in  Region 
lib,  Case  one  we  should  use  the  function  ,4  f[£(y)] t0  ensure  a  bounded  solution  since 

Ai[t,  (y)] - >  0  as  £(y) - >  +°° 

and 

Bi[t,(y)] - >+°°  as  £(y) - >+<*>. 

Following  a  similar  argument  for  Region  TTa,  Case  two  we  should  only  use  the  function 
Ai[t,(y)].  Thus,  by  discarding  the  Bairy  function  in  Region  lib,  Case  one  and  Region  Ila, 
Case  two  we  have  ensured  that  our  solution  will  remain  bounded  above  and  below  the 
source  as  £(y)  approaches  positive  infinity. 

Second,  for  Region  Ila,  Case  one  and  Region  lib,  Case  two  the  value  of  n2(y)  is 
greater  than  one.  Thus  k2(y )  will  take  on  positive  values,  which  leads  to  negative  values 
for  £(y).  Therefore,  both  j4/[£(y)]  and  2?/[£(y)]  should  be  used  in  Region  Ila,  Case  one 
and  Region  lib,  Case  two  since  both  the  Airy  and  Bairy  functions  are  well  behaved  for 
negative  arguments.  The  asymptotic  forms  of  the  Airy  and  Bairy  functions  are  defined  as 
follows,  [Ref.  4,  Chapter  7] 

^i[CCv)]  =  -^[-C(y)]Tsiti[|[-C(y)J5+^]  -Coo  — > +°°.  (2.47) 

and 
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fil‘[C(y)]  =  -^r[-CCv) 4  c°s[f[-C(y)]J  +  J]  ~ ^(y)  — >  +°°.  (2.48) 


Thus,  to  determine  the  general  solution  of  the  inhomogeneous  depth  equation  we 
must  ensure  the  solution  is  bounded  in  Region  lib,  Case  one  and  Region  Ha,  Case  two, 
i.e.,  D  —  0  (Case  one),  see  Eq.(2.42)  and  B  =  0  (Case  two),  see  Eq.(2.41).  Next,  we  must 
ensure  the  solution  is  in  the  form  of  a  propagating  wave  in  Region  Ila,  Case  one  and 
Region  lib,  Case  two,  i.e.,  A  -  By  and  a  =  B/A  (Case  one),  see  Eq.(2.41),  and  C  =  B  and 
b  =  D/C  (Case  two),  see  Eq.(2.42).  In  summary,  we  have  the  following  forms  for  the 
free-space  propagation  solutions  for  each  case  based  on  Eq.(2.41)  and  Eq.(2.42): 


Case  one:  t>rJk,,y)  =  Br  [Ai[t,(y)]  +  aBi[t,(y)]] 

for  y  <  y0 

(2.49) 

O')]] 

for  y  >  y0 

(2.50) 

Case  two:  ■=  Af  Mi[CO)]] 

for  y  <  y0 

(2.51) 

for  y  >  y0 

(2.52) 

where  Ay  and  By  are  arbitrary  constants  and  a  and  b  are  known  constants  to  be  discussed 
later.  Thus,  for  each  case,  we  have  two  equations  and  two  unknowns. 

To  solve  for  the  two  unknown  constants,  we  apply  the  boundary  conditions  of 
continuity  of  acoustic  pressure  and  a  discontinuity  in  the  y  component  of  the  acoustic 
fluid  velocity  vector  at  the  source  depth  y0: 

1 .  continuity  of  acoustic  pressure: 

<bfib(kr,y0)-<i>fla(kr,yo)  =  0  (2.53) 

2.  discontinuity  in  they  component  of  the  acoustic  fluid  velocity  vector: 

^ >/2b(kr,yo )  -  J^<1 >/2a(kr,y0 )  =  d  (2.54) 
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where  d  —  is  the  finite  discontinuity,  [Ref.  2,  Appendix  3C]. 
2jc 

CASE  ONE: 

Substituting  Eq.(2.49)  and  Eq.(2.50)  into  Eq.(2.53)  yields 


AyAi&(yo)]  - By[Ai&(y0)]  +aBi&(y0)]]=  0, 


and  upon  solving  for  Ay  we  obtain 

.  -  tjgCfeaM  +  ggffCCyo)]  3 

y~y  Ai&(y0)} 


(2.55) 


(2.56) 


Substituting  Eq.(2.49)  and  Eq.(2.50)  into  Eq(2.54)  yields 

-ia.i)>AyAi%(yo)]  -  (-)(«,  )^/'[C(yo)]  +  aBi%(y0)]]  =d,  (2.57) 

and  upon  substituting  Eq.(2.56)  into  Eq.(2.57)  and  rearranging  terms,  we  obtain 

-  /«KO'o>Jfl/'[?0’o)]]  = - (2.58) 


Since,  see  [Ref.  1,  Chapter  4] 

BiVi(yo)]Ai%(y o)]  -Ai&(y0)]Bi%(y0)]  =  -£  (2.59) 

which  is  known  as  the  Wronskian  of  the  Airy  and  Bairy  functions,  substituting  Eq.(2.59) 
into  Eq.(2.58)  and  solving  for  By  yields 

By - (2.60) 

a(  ai)j 
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Substituting  Eq.(2.60)  into  Eq.(2.56): 


^  =  — r^/KCyo)]  +  aBi&(yo)]]. 

a(  ai)3 


(2.61) 


:ASE  TWO: 


Substituting  Eq.(2.51)  and  Eq.(2.52)  into  Eq.(2.53)  yields 


By[Ai&(y o)]  +  bBi&(yo)]]  -AyAi&(y0)]  =  0, 


(2.62) 


and  upon  solving  for  A  ,  we  obtain 


A  =B  M/[COo)3  +  ^/[C(vo)]l 

'  '  ^K(yo)l 


(2.63) 


Substituting  Eq.(2.51)  and  Eq.(2.52)  into  Eq.(2.54)  yields 


— <oti)3^M/[C(Vo)]  +  bB^{y0)}}-{-)(a^AyA%(yo)\=d, 


(2.64) 


and  upon  substituting  Eq.(2.63)  into  Eq.(2.64)  and  rearranging  terms,  we  obtain 


Since 


^o)]M/K(yo)3^/[^o)]  "  fi/'fC(yo)^//[COo)]]= - ^ 

(oti)3 


^/TCOo^'lCOo)]  -Bi&(yo)]Ai%<yo)]  =  ±, 


(2.65) 


(2.66) 


substituting  Eq.(2.66)  into  Eq.(2.65)  and  solving  for  By  yields 


By  = - ^-rAi&iyo)]. 

b(aip 


(2.67) 
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Substituting  Eq.(2.67)  into  Eq.(2.63)  yields 


Ay  =  - 


b{ ai)2 


7MCO0)]  +  Wf/Kto)]]. 


(2.68) 


Values  for  the  constants  a  and  b  must  be  determined  before  we  have  a  complete 
solution.  To  determine  the  value  of  a,  recall  Eq.(2.49)  and  Eq.(2.60).  Substituting 
Eq.(2.60)  into  Eq.(2.49)  and  rearranging  terms  yields 


+  frKO)]],  yzy0. 

(ai)3 


(2.69) 


Note  that  (  see  Tables  21  through  23,  located  in  the  Appendix)  as  y  — >  -  00  ,  co)->  —  oo  # 

Recall  the  asymptotic  expansions  of  Aifciy)]  and  #*[CO)]; 


^  =1  ,  3 


^TCO)] = ^KO)l 4  sin[fK(v)]2  +  f]. 


-CO)  ->  +°° 


(2.70) 


£'TCO)]  =  -j=-hC 0)1  4  cos[f  [-£(»]  2  + 


-CO)  — >  +°°- 


(2.71) 


Therefore,  substituting  Eq.(2.70)  and  Eq.(2.71)  into  Eq.(2.69)  yields 


=  —  ,K0'°|)1K(y)1^r^n(fKO')]l  +  f)  +  cosffKWll  +  f)l  (2.72) 
Vrc(oti)3  L  J 


Since  this  derivation  corresponds  to  a  time-harmonic  solution  with  time  dependence 
given  by  exp(+j2itft),  in  order  to  get  <bJ2a(kr  ,y)  into  the  form  of  a  wave  traveling  in  the 
positive  -£(y)  direction,  we  must  choose  a  =  +j  since 


exp(-y'[-CO)])  =  cos([— CO)]) -j  sin  ([-£(»]) . 


(2.73) 
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To  determine  the  value  of  b,  we  follow  a  similar  logic.  Substituting  Eq.(2.67)  into 
Eq.(2.52)  and  rearranging  terms  yields 

+  5/TCO)]],  T-To  •  (2.74) 

(ai)3 

Note  that  (see  Tables  18  through  20,  located  in  the  Appendix)  as  y  — >  +  oo 

Using  the  same  logic  as  before,  since  -£(y)  -->  -H*>,  we  must  pick  b  =  +j  to  ensure  the 
correct  form  of  Eq.(2.74). 

In  summary,  we  have  the  following  formulas: 

CASE  ONE: 

<&/2a(kr,y)  =  nd4l^(-Vo^  [Ai[^(y)  ]  +y5/[^0)]}  ,y<y0  (2.75) 

(+/)(ot  l)3 


and 


<&fib(kr,y) 


nd{Ai\^y 0)]  ^/[C(yo)]} 

(+y')(a,)3 


Ai[£,{y)\  ,y>yo 


(2.76) 


CASE  TWO: 


Q>/2a(kr,y)  =  - 


^M/[C(yo)]  +7^[^(vo)]} 

(-U')(ai)3 


•^/TCCv)]  ,y^yo 


(2.77) 


and 


-  -  —  4 ' ;CV°)1  M/fCCv)]  ,.y>y„ 

(+/)«*,)> 


(2.78) 
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Recall  that  d=-^~. 

2n 

Now  that  a  solution  for  <&f(kr  Vv)  has  been  derived,  all  that  is  left  to  determine  the 
velocity  potential  is  to  insert  Q>f{kr  y)  into  Eq.(2.3)  and  integrate.  It  is  interesting  to 
note,  as  described  in,  [Ref.  2,  Chapter  9],  that: 


*t(*r,y)  =Hm{fJr,y  ;yo )  (2.79) 

where  HJff^y ;  yB )  is  the  time-invariant,  space-variant  ocean  medium  transfer  function. 
Thus,  as  described  in  the  Introduction  of  this  thesis,  Eq.(2.75)  through  Eq.(2.78)  will  be 
added  to  the  already  existing  library  of  ocean  medium  transfer  functions  in  the  program 
LSVOCN.  With  this  derivation  and  new  transfer  function,  we  can  compare  the  results 
from  LSVOCN  to  the  RRA  algorithm  to  determine  the  accuracy  of  the  RRA  algorithm. 
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m.  COMPARISON  OF  THE  "AIRY  FUNCTION"  SOLUTION 
(MAGNITUDE  AND  PHASE) 

WITH  THE  RRA  ALGORITHM’S  SOLUTION 

A.  THE  TIME-HARMONIC  FREE-SPACE  GREEN'S  FUNCTION  AND 
NUMERICAL  ANALYSIS  TECHNIQUES 

1.  The  Time-harmonic  Free-Space  Green's  Function  for  a  Homogeneous 
Ocean  Medium 

Once  the  solution  to  the  inhomogeneous  Helmholtz  equation  for  a  free-space 
propagation  problem  with  the  square  of  the  index  of  refraction  a  linear  function  of  depth 
had  been  derived,  the  solution  was  converted  to  FORTRAN  computer  code  and  added 
to  the  computer  program  LSVOCN  as  a  new  ocean  medium  transfer  function.  LSVOCN 
will  accept  as  primary  inputs: 

1.  Sound-speed  at  the  ocean  surface  and  source,  in  meters/second. 

2.  Receiver  depth  and  horizontal  range  relative  to  the  source,  in  meters. 

3.  Frequency  the  source  will  operate,  in  Hertz. 

The  output  of  LSVOCN  is  either  the  velocity  potential  (magnitude  and  phase)  or 
the  acoustic  pressure  (magnitude  and  phase)  at  a  specified  receiver  location.  To  test  the 
RRA  Algorithm,  the  depth  and  horizontal  range  coordinates  of  selected  points  along  a 
ray  path,  propagated  by  the  RRA  Algorithm,  will  be  input  into  LSVOCN  as  receiver 
coordinates.  The  calculated  acoustic  pressures  (magnitude  and  phase)  from  LSVOCN 
corresponding  to  these  selected  points,  will  be  compared  with  the  calculated  acoustic 
pressures  (magnitude  and  phase)  from  the  RRA  Algorithm. 

Before  beginning  the  testing  process,  we  needed  to  ensure  that  both  the 
theoretical  solution  and  FORTRAN  computer  code  were  correct.  In  order  to  perform  this 
validation,  test  cases  were  developed  for  which  theoretical  values  of  the  velocity 
potential  could  be  calculated  and  compared  with  the  predicted  values  found  using 
LSVOCN.  The  time-harmonic  free-space  Green's  function  which  corresponds  to  a 
spherical  sound  source  operating  in  the  monopole  mode  of  vibration  in  a  homogeneous 
ocean  medium,  was  selected  as  the  basis  for  these  test  cases. 
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As  shown  in  [Ref.  2,  Chapter  4],  the  velocity  potential  of  the  acoustic  field 
produced  by  a  spherical  sound  source  in  the  monopole  mode  of  vibration  operating  in  an 
unbounded,  homogeneous  medium  is  given  by 


<P/<* , e , y)  =  Agyfelg)  ,R>  a. 


where  S0  is  the  source  strength  (volume  flow  rate)  in  mVs, 


is  the  time-harmonic  ffee-space  Green's  function  (exp(jcjt)  factor  suppressed)  and  R  is 

the  spherical  range  between  source  and  receiver.  In  the  limit  as  a  0  (modeling  a  point 
source)  Eq.(3.1)  becomes 


cp/(K,e,V|/)  =  So^rlO)  =  _s0£2LM 


,  R>  0. 


Using  Eq.(3.3),  theoretical  values  for  the  magnitude  and  phase  of  an  acoustic 
field  at  a  receiver  in  an  unbounded,  homogeneous  ocean  medium  located  R  meters  from  a 
point  source  can  easily  be  calculated.  Although  these  values  are  valid  only  for  the 
homogeneous  case,  as  the  modeled  ocean  medium  that  LSVOCN  uses  becomes  more 
homogeneous,  then  the  calculated  velocity  potential  (magnitude  and  phase)  from 
LSVOCN  should  approach  the  predicted  values  calculated  by  Eq.(3.3).  The  theoretical 
magnitude  of  the  velocity  potential  calculated  by  Eq.(3.3)  is  given  by 


M  = 


The  theoretical  phase  of  the  velocity  potential  calculated  by  Eq.(3.3)  is  given  by 
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(3.5) 


Z<p/=Jc(±l  - 

where /is  the  frequency  at  which  the  source  operates,  R  is  the  spherical  range  between 
the  source  and  receiver,  and  c  is  the  speed  of  sound.  In  order  to  test  the  Airy  function 
solution  and  the  corresponding  FORTRAN  computer  code  in  LSVOCN,  Eq.(3.4)  and 
Eq.(3.5)  were  used  in  three  different  Green's  function  test  cases,  defined  as  follows: 

1.  With  the  receiver  and  source  depth  equal,  and  with  the  receiver  one  meter  from 
the  source,  vary  the  speed  of  sound  at  the  source  so  that  it  approaches  the  speed 
of  sound  at  the  surface. 

2.  With  the  receiver  and  source  depth  equal,  and  the  source  sound-speed  constant, 
vary  the  horizontal  range  between  the  source  and  receiver. 

3.  With  the  source  sound-speed  constant,  vary  the  receiver  depth. 

Each  of  these  test  cases  will  be  covered  in  detail  in  Section  B  of  this  chapter. 

2.  Numerical  Analysis  Techniques 

Recall  from  Eq.(2.3)  that 

^>fir,y)  =H~'  {<&f{kr,y)}  =  j~  <j>fikr,y)Jo{krr)krdkr .  (3.6) 

In  order  to  calculate  the  exact  "Airy  function"  solution  at  a  given  receiver  depth  and 
horizontal  range,  the  integral  in  Eq.(3.6)  must  be  computed.  As  we  discovered,  this 
integral  is  difficult  to  evaluate  numerically.  The  integrand  in  Eq.(3.6)  is  oscillatory  due 
to  the  Bessel  function,  when  negative  arguments  of  the  Airy  and  Bairy  functions  are 
involved.  Numerous  overflow  and  underflow  errors  were  generated  when  Eq.(3.6)  was 
calculated  numerically  with  the  computer.  To  eliminate  overflow  errors,  the 
exponentially  scaled  versions  of  the  Airy  and  Bairy  functions  were  used.  The  following 
example  illustrates  the  method  of  using  the  exponentially  scaled  Airy  and  Bairy  functions 
to  eliminate  overflow  errors. 

In  order  to  solve  Eq.(2.75)  or  Eq.(2.78)  for  the  case  of  propagating  waves,  it  is 
necessary  to  multiply  the  Airy  and  Bairy  functions  together.  Although  the  argument  of 
these  functions,  £(y),  is  more  likely  to  take  on  negative  values  in  regions  11a,  Case  one 
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and  region  lib,  Case  two,  once  the  propagating  waves  turn  evanescent,  the  value  for  £(v) 
is  forced  to  take  on  positive  values  which  cause  the  Bairy  function  to  blow  up,  leading  to 
overflow  errors.  By  experimentation  it  was  discovered  that  the  computer  could  handle 
values  of  £(y)  up  to  positive  forty,  without  generating  overflow  errors.  However,  to 
include  most  of  the  evanescent  waves,  £0)  must  be  allowed  to  take  on  much  larger 

positive  values,  which  requires  the  use  of  the  exponentially  scaled  Airy  and  Bairy 
functions: 

1. DAI(x),  DBI(x)  -  double-precision  IMSL  Airy  and  Bairy  functions. 

2. DAIE(x),  DBIE(x)  -  double-precision  IMSL  exponentially  scaled  and  Bairy 
functions. 


where 


DAIE(x)  =  DAI(x )  x  exp(f*2 ) 


DBIE(x )  =  DBl(x)  x  exp(-  \x* ) . 


Therefore  when  multiplying  ^/[^(y)]  and  first  form  the  difference  term: 


£)7/^F=  |[ccv»)  ’  -  CO)  ’  ]- 


Then  perform  the  required  multiplication  as  follows: 


^.[C0o)l  x^fCfy)]  =  exp (-DIFF)  xDAIE&(y0)]  xDBIE[^(y)] .  (3.10) 

This  method  of  multiplying  the  Airy  and  Bairy  functions  together  allows  for  much  larger 
values  of  £(y)  to  be  used.  Thus  most  of  the  evanescent  waves  could  be  included. 
Underflow  errors  did  not  offer  the  same  problems  that  overflow  errors  caused  since  the 
machine  handled  underflow  errors  without  termination  of  the  computer  program 
LSVOCN.  However,  small  inaccuracies  from  underflow  errors  had  to  be  expected. 
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These  same  techniques  can  be  used  to  solve  Eq.(2.76)  and  Eq.(2.77).  Overflow 
errors  were  not  a  problem  for  these  two  equations  since  the  Bairy  function  was  not  part 
of  the  solution.  Graphical  presentations  of  the  Airy  and  Bairy  functions  are  included  in 
Figures  7  and  8.  Tabulated  values  of  the  Airy  and  Bairy  functions,  the  exponentially 
scaled  Airy  and  Bairy  functions,  and  EDAIE(x)  and  EDBIE(x)  are  included  in  Table  1. 
The  functions  EDAIE(x)  and  EDBIE(x)  are  the  exponentially  scaled  Airy  and  Bairy 
functions  multiplied  by  the  appropriate  exponential  factor  that  eliminates  the  exponential 
scaling.  These  values  were  also  tabulated  to  ensure  that  the  numerical  method 
represented  by  Eq.(3.9)  and  Eq.(3.10)  would  produce  the  same  values  as  if  the 
exponentially  scaled  functions  were  not  used.  This  can  be  observed  in  Table  1  by 
noticing  that  DAI(X)  =  EDAIE(X)  and  that  DBI(X)  =  EDBIE(X). 

To  solve  for  the  velocity  potential,  Eq.(3.6)  must  theoretically  be  integrated  from 
zero  to  infinity.  In  practice,  the  upper  limit  of  integration  need  not  be  infinity,  but  rather 
a  sufficiently  large  number  so  that  the  integrand  becomes  negligible  in  value,  in  other 
words,  all  significant  evanesent  waves  are  included.  A  method  was  developed  to 
automatically  generate  the  upper  limit  of  integration.  This  method  will  be  discussed  next. 

We  begin  with  Eq.(2.45),  rewritten  here  for  convenience: 

C(y)  =  -(ai)"3^0).  (3.11) 

Equation  (3.1 1)  can  rewritten  as 

C 00  = - j===[k2(y)  -  k2r],  .  (3.12) 

ko)2 


and  upon  rearranging  terms  and  solving  for  k„  we  obtain 


kr=+Jt2(y)  +  ^(fliA^)2 


coo. 


(3.13) 
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DAIE(X) 


0.0E+00 

0.355028E+00 

0.1E+02 

0.158I24E+00 

0.1E+03 

0.891969E-01 

0.1E+04 

0.501642E-01 

0.1E+05 

0.282095E-01 

0.1E+06 

0.158634E-01 

0.1E+07 

0.892062E-02 

0.1E+08 

0.501643E-02 

0.1E+09 
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0.1E+10 

0.158634E-02 

O.IE+11 

0.892062E-03 
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0.501643E-03 
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0.158634E-03 
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0.1E+16 
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0.892062E-06 

0.1E4  24 

0.501643E-06 

0.1E+25 

0.282095E-O6 
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0.158634E-06 
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0.892062E-07 
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0.501643E-07 

0.1E+29 

0.282095E-07 

0.1E+30 

0.158634E-07 

0.1E+31 

0.892062E-08 

0.1E+32 

0.501643E-08 

0.  lE-*-33 

0.282095E-08 

0.1E+34 

0.158634E-08 

0.1E+35 

0.892062E-09 

0.1E+36 

0.501643E-09 

0.1E+37 

0.282095E-09 

O.IE+38 

0.158634E-09 

0.IE+39 

0.892062E-10 

0.1E+40 

0.501643E-10 

0.1E+41 

0.282095E-10 

0.1EM2 

0.158634E-10 

Airy  and  Exponentially  Scaled  Airy  Functions 


Airy  Function  with  -x  and  +x  Argument  Values 


Figure  7.  Tabulated  values  of  DAIE(x)  and  plots  of  DAI(x),  DAl(-x),  DAIE(x),  and  EDAIE(x). 
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Figure  8.  Tabulated  values  of  DBIE(x)  and  plots  of  DBI(x),  DBI(-x),  DBIE(x),  and  EDBIE(x). 
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X 

DAI(X) 

DAIE(X) 

EDAIE(X) 

I  X 

DBIfX) 

DBEE(X) 

EDBIE(X) 

0.00 

0.355028E+00 

0.355028E+00 

0.355028E+00 

0.00 

0.614927E+00 

0.614927E+00 

0.614927E+00 

0.10 

0.329203E+00 

0.3362 17E400 

0.329203E+00 

0.10 

0.659862E+00 

0.646096E+00 

0.659862E+O0 

0.20 

0.303703E+00 

0.322363E-KK) 

0.303703E+00 

020 

0.705464E+00 

0.66462 8E +00 

0.7054 64E+00 

0.30 

0278806E+00 

0.31 1084E+00 

0278806E-HX) 

0.30 

0.752486E+00 

0.6744 09E+00 

0.752486E+O0 

0.40 

0.254742E-KX) 

0.301541E+00 

0254742E+00 

0.40 

0.801 773E+00 

0.677338E+O0 

0.801 773E+00 

0.50 

0.23 1694 E+00 

0.293277E+O0 

0.23 1694 E+00 

0.50 

0.854277F+O0 

0.674892E+00 

0.854277E+00 

0.60 

0209800E+00 

0286000E-KX) 

0209800E+00 

0.60 

0.911063E+00 

0.668324E+O0 

0.911063E+00 

0.70 

0.189162E+00 

02795 13E+00 

0.189I62E+00 

0.70 

0.973329E+00 

0.658708E+00 

0.973 329E +00 

0.80 

0.169846E+00 

0273670E+00 

0.169846E+00 

0.80 

0.104242E+01 

0.646954E+00 

0.104242E+01 

0.90 

0.151887E+00 

0268364E+00 

0.151887E+00 

0.90 

0.1 1 1987E+01 

0.633817E+00 

0.1 1 1987E+01 

1.00 

0.I35292E+O0 

02635 14E+00 

0.135292E+00 

1.00 

0.120742E+01 

0.619912E+00 

0.120742E+01 

1.10 

0.120049E+00 

0.259052E+00 

0.1 2004 9E+00 

1.10 

0.130707E+OI 

0.605721  E+OO 

0.130707E+01 

1.20 

0.106126E+00 

0254928E+00 

0.106126E+00 

120 

0.142113E+01 

0.591 614E+00 

0.1421 13E+01 

1.30 

0.934746E-01 

0251098E+00 

0.934746E-01 

1.30 

0.155228E+01 

0.577859E+00 

0.155228E+01 

1.40 

0.8203  80E-01 

0247527E+00 

0.8203  80E-0I 

1.40 

0.170366E+01 

0.564 646E+O0 

0.170366E+01 

1.50 

0.7I7494E-01 

0244185E+00 

0.717495E-01 

1.50 

0.187894E+01 

0352094E+00 

0.187894E+01 

1.60 

0.625369E-01 

0241048E-KK) 

0.625369E-01 

1.60 

0208247E+01 

0.540272E+00 

0208247E+01 

1.70 

0.543248E-01 

0238094E+00 

0.543248E-01 

1.70 

0231941E+01 

0.529208E+00 

0231941E+01 

1.80 

0.470362E-01 

0235306E+00 

0.470362E-01 

1.80 

0259587E+01 

0.518898E+00 

0259587E+01 

1.90 

0.405944E-01 

0232667E+00 

0.405944E-01 

1.90 

0291918E+01 

0309320E+00 

0291918E+01 

2.00 

0.34924  IE-01 

0230165E+00 

0.34924  IE-01 

2.00 

0.329809E+01 

0.500437E+00 

0.329809E+01 

2.10 

0299526E-01 

0.227786E+O0 

0299526E-01 

2.10 

0.374315E+01 

0.492203E+00 

03743 15E+01 

220 

0256104E-01 

0225522E+00 

0256104E-01 

220 

0.426704E+01 

0.484567E+00 

0.426704E+01 

2.30 

0218320E-01 

0223361E-+00 

0218320E-01 

2.30 

0.488506E+01 

0.477480E+00 

0.488506E+01 

2.40 

0.185561E-01 

0221297E+00 

0.185561E-01 

2.40 

0.561 577E+01 

0.470890E+00 

0.561577E+01 

2.50 

0.157259E-01 

0219322E+00 

0.157259E-01 

2.50 

0.648166E+01 

0.464750E+00 

0.648166E+01 

2.60 

0.132893E-01 

02I7429E+00 

0.132893E-01 

2.60 

0.751009E+01 

0.45901 7E +00 

0.751009E+01 

2.70 

0.1 1 1985E-01 

0.2I5613E+00 

0.1 1 1985E-01 

2.70 

0.873439E+OI 

0.453648E+00 

0.873439E+01 

2.80 

0.941050E-02 

0213868E+00 

0.941 05  IE-02 

2.80 

0.101953E+02 

0.44 8 60 8E +00 

0.101953E+02 

2.90 

0.78863  IE-02 

02121 89E400 

0.78863  IE-02 

2.90 

0.1 19426E+02 

0.443863E+00 

0.1 19426E+02 

3.00 

0.6591 14E-02 

0210572E+00 

0.6591 14E-02 

3.00 

0.140373E+02 

0.439384E+00 

0.140373E+02 

3.10 

0.549399E-02 

020901 3E+O0 

0.549399E-02 

3.10 

0.165547E+02 

0.435 146E +00 

0.165547E+02 

320 

0.456744E-02 

0207509E+00 

0.456744E-02 

320 

0.195870E+02 

0.431 125E+00 

0.195870E+02 

3.30 

0.378729E-02 

0206056E+00 

0.378729E-02 

3.30 

0232483E+02 

0.427301E+00 

0J232483E+O2 

3.40 

0.313234E-02 

020465 1E+00 

0.313234E-02 

3.40 

0276796E+02 

0.423657E+00 

0276796E+02 

3.50 

02584 10E-02 

0203292E+00 

0258410E*02 

3.50 

0.330555E+02 

0.4201 77E +00 

0.330555E+O2 

3.60 

0212648E-02 

0201976E4O0 

0212648E-02 

3.60 

0.3 95 92 7E +02 

0.416848E+00 

0.395927E+02 

3.70 

0.1 74557E-02 

0.200700E+00 

0.I74557E-02 

3.70 

0.475607E+02 

0.41 365 6E  +-00 

0.475607E+02 

3.80 

0.1 4293  9E-02 

0.199462E+00 

0.142939E-02 

3.80 

0.572954E+02 

0.410592E+00 

0.572954E+O2 

3.90 

0.1 16765E-02 

0.198261E+00 

0.1 16765E-02 

3.90 

0.692 160E +02 

0.407646E+00 

0.6921 60E +02 

4.00 

0.95 1564E-03 

0.197095E+00 

0.951564E-03 

4.00 

0.838471E+O2 

0.404  809E +00 

0.83  8471E+02 

4.10 

0.773630E-03 

0.I95961E+OO 

0.773629E-03 

4.10 

0.101 846E+03 

0.402075E+00 

0.101846E+03 

420 

0.627496E-03 

0.194859E+00 

0.627496E-03 

420 

0.124038E+03 

0.399435E+00 

0.124038E+03 

4.30 

0.507787E-03 

0.193786E+00 

0.507787E-03 

4.30 

0.151462E+O3 

0.3968 84E+00 

0.151462E+03 

4.40 

0.409974E-03 

0.1 9274 1E+O0 

0.409974E-03 

4.40 

0.185428E+03 

0.3944 17E +00 

0.1 8542  8E +03 

4.50 

0.330250E-03 

0.191 724E-KX) 

0.3  302 5 OE -03 

4.50 

0227588E+03 

0.3 9202 7E +00 

0227588E+03 

4.60 

0265432E-03 

0.190732E+00 

0265432E-03 

4.60 

02  8003 6E +03 

0.389712E+00 

0.280036E+03 

4.70 

0212861E-03 

0.189765E+00 

0212861E-03 

4.70 

0.345426E+03 

0.387466E+00 

0.345426E+03 

4.80 

0. 1 70326E-03 

0.188822E+00 

0.170325E-03 

4.80 

0.427 126E +03 

0.385286E+00 

0.427126E+O3 

4.90 

0.135992E-03 

0.187901E+00 

0.135992E-03 

4.90 

0.529425E+03 

0.383168E+00 

0.529425E+03 

5.00 

0.108344E-03 

0.187002E+00 

0.108344E-03 

5.00 

0.657792E+03 

0381 109E+00 

0.657792E+03 

5.10 

0.861 324E-04 

0.186124E+00 

0.861 324E-04 

5.10 

0.819209E+03 

0.379105E+00 

0.8192J0E+03 

520 

0.683285E-04 

0.185265E+00 

0.683285F.-04 

520 

0.1 02261  E+<M 

0377155E+00 

0.102262E+04 

5.30 

0.540905E-04 

0. 1 84426E+00 

0.540905E-04 

5.30 

0.12794 7E +04 

0.375256E+00 

0.127947E+04 

5.40 

0.427299E-04 

0.183605E+00 

0.427298E-04 

5.40 

0.160448E+04 

0.373405E+00 

0.160448E+04 

5.50 

0.336853E-04 

0.182802E+00 

0.336853E-04 

5.50 

0201658E+04 

0.371600E+00 

0.201658E+04 

5.60 

0265006H-04 

0.182015E+00 

0.265006E-04 

5.60 

0254018E+04 

0.369839E+00 

0.25401 8E+04 

5.70 

0208058E-04 

0.181246E+00 

020805  8E-04 

5.70 

0.320680E+04 

0.368120E+00 

0.320680F.+04 

5.80 

0.16301 7E-04 

0.1 80491 F-KK) 

0.16301 7E-04 

5.80 

0.405 720E +04 

0366441 E+00 

0.405 720E +04 

5.90 

0.127471  E-04 

0.179753E+00 

0.127471E-04 

5.90 

0.5  J4421E  >04 

0.364800E  -*00 

0.5 14422E  >04 

6.00 

0.994769E-05 

0.179028E+00 

0.994  769E-05 

6  00 

0.653645E+04 

0.363197E+00 

0.653645E+O4 

6.10 

0.774773E-05 

0.1783 1 8E+00 

0.774773E-05 

6.10 

0.832309E+04 

0.361629E*00 

0.832309EKM 

620 

0.602246E-05 

0.177622E+00 

0.602246E-05 

620 

0. 1 06204  E +05 

0.360095E+00 

0.106204F.+05 

6.30 

0.467226E-05 

0.176939E+00 

0.467226E-05 

630 

0.135799E+O5 

0358593E+00 

0.1 35800b +05 

Table  1.  Tabulated  values  for  DAI(x),  DAIE(x),  EDAIE(x),  DBl(x),  DBIE(x),  and  EDBIE(x). 
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Before  the  exponentially  scaled  Airy  and  Bairy  functions  were  used,  the  largest 
positive  value  of  the  argument  for  the  Bairy  function  that  could  be  evaluated  by  the 
computer  without  causing  overflow  errors  was  £(y)  =  =  40.  By  using  the 

exponentially  scaled  Airy  and  Bairy  functions,  it  was  discovered  that  could  take  on 
extremely  large  values,  for  example,  =  1012  -  1016.  Substituting  ^  into  Eq.(3.13) 
for  £(y)  resulted  in  the  following  equation  for  the  upper  limit  of  integration. 


*rmaX=+^2(y)  +  £ 


(3.14) 


Therefore,  the  upper  limit  of  integration  used  in  Eq.(3.6)  could  automatically  be  adjusted 
to  the  maximum  possible  value  regardless  of  the  frequency  selected.  However,  due  to  the 
oscillatory  nature  of  the  Bessel  function  and  the  dependence  of  its  argument  on  kr  and  r, 
the  integrand  of  Eq.(3.6)  becomes  increasingly  difficult  to  integrate  numerically  as  kr 
increases  (for  a  fixed  r),  because  as  the  wave  number  of  the  Bessel  function  increases,  the 
function  becomes  more  highly  oscillatory.  To  ensure  accurate  results  from  the  computer 
code,  the  region  of  integration  had  to  be  divided  into  subintervals  in  performing  the 
numerical  integration.  This  method  enabled  LSVOCN  to  calculate  increasingly  many 
evanescent  waves  providing  more  accurate  results. 

B.  COMPARISON  OF  LSVOCN  WITH  THE  TIME-HARMONIC 
FREE-SPACE  GREEN’S  FUNCTION  FOR  A  HOMOGENOUS  OCEAN 
MEDIUM 

1.  Green’s  Function  Test  Case  One  -  Vary  the  Source  Sound-Speed 

The  first  test  case  involved  varying  the  source,  operating  at  frequencies  of 
1000  Hz,  250  Hz,  and  50  Hz,  sound-speed  while  maintaining  the  surface  sound-speed 
fixed,  using  both  positive  and  negative  gradients  for  the  square  of  the  index  of  refraction. 
Each  case  had  an  initial  source  sound-speed  of  1515  m/s  for  the  negative  gradient  of  the 
square  of  the  index  of  refraction,  and  1485  m/s  for  the  positive  gradient  case.  The  source 
and  receiver  depths  were  held  constant  at  150  meters,  with  the  receiver  located  one  meter 
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in  horizontal  range  from  the  source.  For  each  successive  run  of  the  test  case,  the  source 
sound-speed  was  decreased/increased  by  a  factor  of  3  m/s  until  1501  m/s  and  1499  m/s 
was  reached,  respectively,  while  maintaining  the  surface  sound-speed  constant  at  1500 
m/s.  Thus,  the  source  sound-speed  was  never  allowed  to  take  on  a  value  of  1500  m/s, 
ensuring  that  the  Airy  function  solution  of  LSVOCN,  which  assumed  a  sound-speed 
profile  where  the  square  of  the  index  of  refraction  was  a  linear  function  of  depth,  would 
remain  valid  while  at  the  same  time  forcing  the  modeled  ocean  medium  to  become  more 
homogeneous.  If  the  derived  Airy  function  solution  and  the  corresponding  computer  code 
in  LSVOCN  were  correct,  then  as  the  source  sound-speed  approached  the  surface 
sound-speed,  the  velocity  potential  (magnitude  and  phase)  values  calculated  by  LSVOCN 
would  converge  approximately  to  the  velocity  potential  (magnitude  and  phase)  values 
predicted  by  Eq.(3.4)  and  (3.5),  respectively.  Because  we  were  making  calculations  at  the 
source  depth,  the  DIFF  term  given  by  Eq.(3.9)  was  equal  to  zero,  which  made  the 
exp (-DIFF)  term  in  Eq.(3.10)  equal  to  one.  This  condition  caused  the  oscillations  of  the 
integrand  in  Eq.(3.6)  to  decay  very  slowly,  thus  requiring  longer  integration  causing  an 
increase  in  CPU  time.  A  summary  of  these  test  cases  is  included  as  Tables  2  through  7. 

The  Delta's  and  percentage  differences  included  in  Tables  2  through  7  were 
calculated  as  follows: 

1.  Delta  magnitude  =  LSVOCN  magnitude  -  Green’s  function  magnitude 

2.  Delta  phase  =  LSVOCN  phase  -  Green's  function  phase 


3.  %  Difference  magnitude  = 


LSVOCN  mag  -  Greoi’s  function  mag 
Green’s  function  mag 


x  100% 


4.  %  Difference  phase  = 


LSVOCN  phase  -  Green’s  function  phase 
Green’s  function  phase 


x  100% 


The  data  in  Tables  2  through  7,  show  that  as  the  source  sound-speed  approaches 
the  surface  sound-speed,  making  the  modeled  ocean  medium  more  homogeneous,  the 
calculated  velocity  potential  (magnitude  and  phase)  values  from  LSVOCN  approach  the 
velocity  potential  (magnitude  and  phase)  values  predicted  by  the  Green’s  function  -  our 
first  validation. 
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|  SOURCE  FREQUENCY  =  1  QOOHz  1 

I  SOURCE  DEPTH  =  RECEIVER  DEPTH  =  1 50  METERS  1 

1  SURFACE  SOUND-SPEED  =  1 500  m/s  | 

Sound  Speed 
Source 

m/s 

Magnitude 
Green’s  Func 

m  2/s 

Phase 

Green’s  Func 
Deg. 

Magnitude 

LSVOCN 

m Vs 

Phase 

LSVOCN 

Deg. 

Deha 

Magnitude 

m2/s 

Delta 

Phase 

Deg. 

%  Difference 
Magnitude 

%  Difference 
Phase 

1515 

0.0795775 

302.37624 

0.0787327 

302.557 

-0.0008448 

0.18076 

-1.0616066 

0.0597798 

1512 

0.0795775 

301.90476 

0.0787842 

302.059 

-0.0007933 

0.15424 

-0.9968898 

0.051089 

1509 

0.0795775 

301.43141 

0.0788504 

301.559 

-0.0007271 

0.12759 

-0.9137005 

0.042328 

1506 

0.0795775 

300.95618 

0.0789223 

301.039 

-0.0006552 

0.08282 

-0.8233483 

0.027519 

1503 

0.0795775 

300.47904 

0.0790248 

300.502 

-0.0005527 

0.02296 

-0.6945431 

0.0076411 

Table  2.  Comparison  of  Green's  Function  with  LSVOCN  for  Test  Case  1 
(Vary  the  Source  Sound-Speed)  -  n2(y)  with  Negative  Gradient. 


I  SOURCE  FREQUENCY  =  250Hz 

|  SOURCE  DEPTH  =  RECEIVER  DEPTH  =  1 50  METERS 

1  SURFACE  SOUND-SPEED  =1500  m/s 

Sound  Speed 
Source 

m/s 

Magnitude 
Green's  Func 

m2/s 

Phase 

Green's  Func 
Deg. 

Magnitude 

LSVOCN 

m2/s 

Phase 

LSVOCN 

Deg. 

Delta 

Magnitude 

mVs 

Deha 

Phase 

Deg. 

%  Difference 
Magnitude 

%  Difference 
Phase 

1515 

|  0.0795775  |  120.59406  |  0.0791376  |  120.992  |  -0.0004399  |  0.39794  |  -0.5527944  |  0.3299831  | 

1512 

0.0795775 

120.47619 

0.0791609 

120.83 

-0.0004166 

0.35381 

-0.5235148 

0.2936763  I 

1509 

0.0795775 

120.35785 

0.0792156 

120.704 

-0.0003619 

0.34615 

-0.4547768 

0.2876007  I 

1506 

0.0795775 

120.23904 

0.07923 

120.502 

-0.0003475 

0.26296 

-0.4366812 

0.2186977  1 

1503 

0.079577 5 

MESEZM 

0.0793777 

120.422 

-0.0001998 

0.30224 

-0.251076 

ammm 

1501 

0.0795775 

120.03997 

0.0794475 

120.263 

-0.00013 

0.22303 

-0.1633628 

0.1857964  ■ 

Table  3.  Comparison  of  Green’s  Function  with  LSVOCN  for  Test  Case  1 
(Vary  the  Source  Sound-Speed)  -  n2(y)  with  Negative  Gradient. 


■  SOURCE  FREQUENCY  =  50Hz  1 

I  SOURCE  DEPTH  =  RECEIVER  DEPTH  =  1 50  METERS  1 

1  SURFACE  SOUND-SPEED  =  1 500  m/s  1 

Sound  Speed 
Source 

m/s 

Magnitude 
Green's  Func 

m2/s 

Phase 

Green's  Func 
Deg. 

Magnitude 

LSVOCN 

m2/* 

Phase 

LSVOCN 

Deg. 

Delta 

Magnitude 

m2/s 

Deha 

Phase 

Deg. 

•/•Difference 

Magnitude 

%  Difference 
Phase 

1515 

0.0795775 

168.11881 

0.079782 

168.366 

0.0002045 

0.24719 

0.2569822 

0.1470329 

1512 

0.0795775 

168.09524 

0.0797383 

168.322 

0.0001608 

0.22676 

0.2020672 

0.1348997 

1509 

0.0795775 

168.07157 

e msm 

168.286 

0.000205 

0.21443 

0.2576105 

0.1275826 

1506 

0.0795775 

168.04781 

0.079723 

168.231 

0.0001455 

0.18319 

0.1828406 

0.1090106 

1503 

0.0795775 

168.02395 

0.0797083 

168.172 

0.0001308 

0.14805 

0.1643681 

0.0881124 

^ \ 

■itiVZj!tV>U 

Table  4.  Comparison  of  Green's  Function  with  LSVOCN  for  Test  Case  1 
(Vary  the  Source  Sound-Speed)  -  n2(y)  with  Negative  Gradient. 
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Sound  Speed  Magnitude  Phase 

Source  Green’s  Func  Green’s  Func 
m/s  mVs  Deg. 


SOURCE  FREQUENCY  =  1000Hz 
SOURCE  DEPTH  =  RECEIVER  DEPTH  =  1 50  METERS 
SURFACE  SOUND-SPEED  =  1500  m/s 


Magnitude 

LSVOCN 

mVs 


Phase 

LSVOCN 

Deg. 


Delta 

Magnitude 

mVs 


0.0795775  I  297.57576  I  0.0787358 


0.0795775  I  298.06452  |  0.0787811 


297.709  |  -0.0008417 


%  Difference  I  %  Difference 
Magnitude  I  Phase 


-1.057711 


0.0447752 


298.174  I  -0.0007964  I  0.10948  I  -1.0007854  1  0.0367303 


0.0795775  j  298.55131  j  0.0788484  I  298.647  I  -0.0007291  I  0,09569  I  -0  9162138  _l  0.0320514 

0.0795775  i  299.03614  j  0.0789214  I  299.099  I  -0.0006561  I  0.06286  I  -0.8244793 


0.0795775  I  299.51904  I  0.079039 


0.0795775  I  299.83989  I  0.079143 


0.0210209 


299.546  1  -0.0005385  1  0.02696  1-0.6766988  1  0.0090011 


299.839  |  -0.0004345  I  -0.00089  |  -0.5460086  I  -0.0002968 


Table  5.  Comparison  of  Green’s  Function  with  LSVOCN  for  Test  Case  1 
(Vary  the  Source  Sound-Speed)  -  n2(y)  with  Positive  Gradient. 


Sound  Speed 
Source 
m/s 


SOURCE  FREQUENCY  =  250Hz 
SOURCE  DEPTH  =  RECEIVER  DEPTH  =  150  METERS 
SURFACE  SOUND-SPEED  =  1 500  m/s 


Magnitude  Phase  Magnitude 

Green's  Func  Green's  Func  LSVOCN 
m’/s  Deg.  ■>'/• 


Phase 

LSVOCN 

I*g 


0.0795775 


0.0795775 


0.0795775 


0.0795775 


0.0795775 


0.0795775 


119.39394  I  0.0791194 


119.51613 


0.0791708 


119.63783 


0.0791928 


119.75904  I  0.0792593 


119.87976  I  0.0793556 


119.95997  I  0.0794502 


Delta 

Magnitude 

mVs 


-0.0004581 


-0.0004067 


-0.0003847 


-0.0003182 


-0.0002219 


— wai 


0.22503  I  -0.1599698 


%  Difference  I  */*  Difference 
Magnitude  I  Phase 


0.3183244 


0.3086362 


0.267616 


0.2513046 


0.221255 


0.1875876 


CMldJI 


Table  6.  Comparison  of  Green's  Function  with  LSVOCN  for  Test  Case  1 
(Vary  the  Source  Sound-Speed)  -  n2{y)  with  Positive  Gradient. 


SOURCE  FREQUENCY  =  50Hz 
SOURCE  DEPTH  =  RECEIVER  DEPTH  =  1 50  METERS 
_ _ SURFACE  SOUND-SPEED  =  1500  m/s 


Magnitude  Phase  Magnitude  Phase 

Green’s  Func  Green's  Func  LSV  OCN  LSVOCN 

mVs  Deg. 


Deha 

Magnitude 

m2/s 


0.0795775 


0.0795775 


0.0795775 


0.0795775 


0.0795775 


0.0795775 


167.87879 


167.90323 


167.92757 


167.95181 


167.97595 


167.99199 


0.0797431 


0.0797786 


0.0797651 


0.0797772 


0.0797195 


168.123 


168.138 


168  136 


168  142 


168  135 


0.0001656 


0.0002011 


0.0001547 


0.0001876 


0.0001997 


0.000142 


0.15905 


0.11101 


%  Difference  j  %  Difference 
Magnitude  I  Phase 


0.145468 


0.1398246 


0.124119 


0.1132408 


0.0946862 


0.0660805 


0.235745 


0.2509503 


0.1784424 


Table  7.  Comparison  of  Green’s  Function  with  LSVOCN  for  Test  Case  1 
(Vary  the  Source  Sound-Speed)  -  n2(y)  with  Positive  Gradient. 


2.  Green’s  Function  Test  Case  Two  -  Vary  the  Horizontal  Range 

The  second  test  case  involved  varying  the  horizontal  range  while  maintaining  the 
source  and  receiver  depths  fixed  at  150  meters.  For  this  case  the  surface  and  source 
sound-speeds  are  not  varied.  If  the  receiver  is  fixed  at  the  source  depth  (150  meters)  and 
the  horizontal  range  is  increased,  then  the  magnitude  of  the  velocity  potential  should 
decrease  by  a  factor  of  l/4jtfl,  independent  of  frequency,  as  predicted  by  Eq.(3.4). 
Furthermore,  it  should  be  observed  that  the  phase  varies  as  predicted  by  Eq.(3.5): 

Z<p/=JC^±1  -  (3.15) 

The  above  observations  are  reasonable  assumptions  to  make  even  though  the  modeled 
ocean  medium  is  not  homogeneous.  As  long  as  the  receiver  is  maintained  at  the  same 
depth  as  the  source,  then  the  sound  speed  is  constant,  which  simulates  a  homogeneous 
ocean  medium. 

To  conduct  the  second  Green's  function  test  case,  the  receiver  depth  and  source 
depth  are  both  fixed  at  150  meters.  The  surface  sound-speed  is  fixed  at  1500  m/s  and  the 
source  sound-speed  is  fixed  at  1515  m/s  for  a  negative  gradient  for  the  square  of  the 
index  of  refraction,  and  1485  m/s  for  a  positive  gradient  for  the  square  of  the  index  of 
refraction.  Horizontal  range  values  selected  were  2.0,  3.0,  and  5.0  meters.  These  values 
were  selected  based  on  the  amount  of  CPU  time  required  to  calculate  the  velocity 
potential  (magnitude  and  phase)  for  each  selected  receiver  depth,  horizontal  range  pair. 
Because  we  were  making  calculations  at  the  source  depth,  the  DIFF  term  given  by 
Eq.(3.9)  was  equal  to  zero,  which  made  the  exp  {-DIFF)  term  in  Eq.(3.10)  equal  to  one. 
This  condition  caused  the  oscillations  of  the  integrand  in  Eq.(3.6)  to  decay  very  slowly, 
thus  making  integration  more  difficult  resulting  in  more  CPU  time  being  required.  In 
addition,  as  kr  increases  (for  a  fixed  r)  the  zeroth-order  Bessel  function  in  Eq.(3.6) 
oscillates  more  rapidly,  making  the  numerical  integration  more  difficult.  Any  calculation 
requiring  more  than  60  minutes  of  CPU  time  was  stopped  by  the  computer.  The  test  case 
results  are  summarized  in  Tables  8  and  9.  Calculations  that  required  greater  than  60 
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minutes  of  CPU  time  are  denoted  by:  >60  min.  The  delta's  and  percentage  differences 
were  calculated  in  the  same  manner  as  discussed  previously.  The  data  in  Tables  8  and  9 
show  that  as  the  horizontal  range  increases  while  maintaining  the  receiver  at  the  source 
depth,  the  velocity  potential  (magnitude  and  phase)  calculated  by  LSVOCN  behaves 
according  to  the  Green's  function  for  a  ffee-space  homogeneous  ocean  medium  as 
predicted  by  Eq.(3.4)  and  Eq.(3.5)  with  S0  =  1  -  our  second  validation. 


SOURCE  DEPTH  =  RECEIVER  DEPTH  =  150  METERS 
SURFACE  SOUND-SPEED  =  1500  m/s  SOURCE  SOUND-SPEED  = 


Delta 

Magnitude 

mV* 


1515  m/s 
% 

Difference 

Magnitude 


Green's  Green's  LSVOCN  LSVOCN 
Func  Fane  Magnitude  Phase 

Magnitude  Phase  mVs  Deg. 

mVs  Deg. 


0.0397887 1 64.752475 1  >60  min  I  >60  min 


0.0265258)  187.12871 1  >60  min  I  >60  min 


0.0 159 155 1 7 1.881 188 1  >60  min  I  >60  min 


0.0397887 161.1881 19  |o.03959l9 1  61.221  |-0.0001968  |  0.032881 1 -0.4946128 


0.0265258 1  1.7821 782 1 0.0264635 


0.0159155)  242.9703  |  >60  min  |  >60  min 


0.0397887|  156.23762  |o.0397634 


%  Difference 
Phase 


0.0537376 


-0.0000623|  0.56082 1 8|  -0Z348657|  31.4683346 


-0.0000253|  0.43238| -0.0635859|  0.2767451 


0.026S258|l44.35644|0.0264378|  144.937  |  -O.OOOQ88|  0.580561  -0.3317525 


0.01591551 120.594Q6|o.0157484i  121.152  I -0.000167ll  0.55794 


Table  8.  Comparison  of  Green's  Function  with  LSVOCN  for  Test  Case 
(Vary  Horizontal  Range )  -  n2(y)  with  Negative  Gradient. 


0.4021712 


0.4626596 


SOURCE  DEPTH 
SURFACE  SOUND-SPEED  _ 

Frequency  Range  Green's  Green's 

Hertz  meters  Fnnc  Ftmc 

Magnitude  Phase 

mVs  Deg. 


0.0397887  I  55.151515 


0.026522 6  |  172.72727 


=  RECEIVER  DEPTH  =  150  METERS 
-  1500  m/s  SOURCE  SOUND-SPEED  -  1485  m/s 


LSVOCN  LSVOCN  Hell* 

Magnitude  Phase  Magnitude 

ml/s  Deg.  m2/s 


>60  min 


>60  min 


Delta  % 

Phase  Difference 

Deg.  Magnitude 


% 

Difference 

Phase 


0.0159155  |  47.878788  |  >60  min 


0.0397887  |  58.787879  )  0.0396074  )  58.825  |-O.OOQ1813|  0.037 1 2 1 1  -0.455657 


0.0265226  |  358.18182  |  0.0264458  |  358.687  |  -0.0000768 1  0.505181-0.2895644 


0.0159155  I  236.9697  |  >60  min  |  >60  min  |  N/A  |  N/A  )  N/A 


0.0397887  I  155.75758  I  0.0397893  1  ^156.212  |  0.0000006 1  Q. 45442)  0.001508 


0.0265226  |  143.63636  |  0.0264342  |  144.218  |-O.OOOQ884|  0.58 164|  -0.3333007 


0.0159155  I  119.39394  I  0.0157441  | 1 19.928  I-00QO1714)  0.534061-1.0769376 


Table  9.  Comparison  of  Green's  Function  with  LSVOCN  for  Test  Case  2 
(Vary  Horizontal  Range  )  -  w2(y)  with  Positive  Gradient. 


0.063144 


0.1410401 


0.2917482 


0.4049393 


0.4473091 
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As  a  final  note  on  the  second  Green's  function  test  case,  as  can  be  observed  from 
Tables  8  and  9,  the  delta  magnitudes  all  agree  out  to  three  significant  digits,  and  in  some 
cases  better.  This  fact  was  consistent  with  the  magnitude  differences  as  calculated  in  the 
first  Green's  function  test  case.  Also,  for  small  velocity  potential  phases,  small  phase 
errors  result  in  large  percentage  differences,  as  expected. 

3.  Green’s  Function  Test  Case  Three  -  Vary  the  Receiver  Depth 

The  third  Green's  function  test  case  involved  varying  the  receiver  depth  above 
and  below  the  source  depth,  while  maintaining  the  source  and  surface  sound-speeds  fixed 
at  1515/1485  m/s  and  1500  m/s,  respectively.  As  was  stated  before,  the  modeled  ocean 
medium  is  not  homogeneous;  however,  for  the  previous  two  Green's  function  tests,  the 
receiver  depth  was  maintained  at  the  source  depth  and  the  source  sound-speed  was  varied 
to  more  closely  match  the  surface  sound-speed  (Test  Case  1)  or  maintained  constant 
(Test  Case  2)  which  simulated  a  homogeneous  ocean  medium.  For  this  case,  the  receiver 
was  moved  slightly  above  and  below  the  source  depth,  and  as  a  result,  the  Green's 
function  results  are  less  valid.  As  can  be  seen  from  Figure  9,  even  though  the  receiver  is 
not  located  at  the  source  depth,  the  actual  change  in  the  sound  speed  at  the  receiver  is 
very  small.  Because  the  ocean  medium  acts  very  much  like  a  homogeneous  medium 
when  receivers  and  sources  are  "close",  the  velocity  potential  (magnitude  and  phase) 
computed  using  the  free-space  Green’s  function  for  a  homogeneous  ocean  medium  is  a 
valid  approximation.  This  test  case  allowed  us  to  validate  the  remaining  part  of  the 
computer  code  in  LSVOCN  for  a  source  and  receiver  at  different  depths,  prior  to 
comparing  it  with  the  RRA  Algorithm.  This  test  case  was  conducted  with  a  source 
operating  at  1000  Hz  (excluding  cases  for  the  source  operating  at  250  Hz  and  50  Hz  since 
validation  at  one  frequency  was  adequate).  The  third  Green's  function  test  case  is 
summarized  in  Table  10.  The  test  case  was  conducted  four  times: 

1 .  For  the  square  of  the  index  of  refraction  with  a  negative  gradient  and  a 
receiver  depth  of  150.5  meters  and  a  horizontal  range  of  0.8660254  meters. 

2.  For  the  square  of  the  index  of  refraction  with  a  negative  gradient  and  a  receiver 
depth  of  149.5  meters  and  a  horizontal  range  of  0.8660254  meters. 


37 


3.  For  the  square  of  the  index  of  refraction  with  a  positive  gradient  and  a  receiver 
depth  of  150.5  meters  and  a  horizontal  range  of  0.8660254  meters. 

4.  For  the  square  of  the  index  of  refraction  with  a  positive  gradient  and  a  receiver 

depth  of  149.5  meters  and  a  horizontal  range  of  0.8660254  meters. 


Region  II:  Receiver  below  the  source 


Figure  9.  Geometry  for  Green's  Function  Test  Case  3  (Vary  Receiver  Depth). 


The  surface  sound-speed  is  maintained  at  1500  m/s.  The  source  sound-speed  is 
maintained  at  1515  m/s  when  the  square  of  the  index  of  refraction  has  a  negative 
gradient,  and  1485  m/s  when  the  square  of  the  index  of  refraction  has  a  positive  gradient. 
Calculations  are  based  on  Eq.(3.4)  and  Eq.(3.5)  where  R  and  c  are  treated  as  variables. 
As  can  be  observed  from  Table  10,  the  change  in  the  sound  speed  at  the  location  of  the 
receiver  is  very  small,  making  our  assumptions  for  this  test  case  reasonable.  Analysis  of 
the  data  presented  in  Table  10  show  that  LSVOCN  calculated  values  for  the  velocity 
potential  (magnitude  and  phase)  are  approximately  equal  to  those  predicted  by  the 
Green's  function  -  our  third  validation.  The  Delta's  and  percentage  differences  were 
calculated  in  the  same  manner  as  discussed  before. 


SURFACE  SOUND-SPEED  =  1500  M/S 

_ _ _  r  =  0.8660254  METERS,  FREQ  =  1000  Hz _ 

Source  Receiver  Sound  Receiver  Green's  Green's  1SVOCN  LSVOCN  Deha  Delta  %  Diff  %DilT 

Sound  Speed  Depth  Func  Func  Mag  Phase  Mag  Phase  Mag  Phase 

Speed  mis  meters  Magnitude  Phase  m:/s  deg.  m’/s  deg. 

nr/s  mVs  deg. 

ISIS  1514.9492501  149.5  0,0795775  302.3683  0.0795773  302.37  -0,0000002  0,0017  -0,0002513  0.0005622 

J51^  J5110507551.  ^150^  00795775  302.3842  0.079575  302.381  -0.0000025  -0.0032  -0.0031416  -0.0010583 

1485  ^485^0492549  J49^  0.0795775  297.5838  0.079575  297.58  -0.0000025  -0,0038  -0.0031416  -0.001277 

1485  |  1484.95075  |  1505  0.0795775  297.5677  0,079577  297.569  -0,0000005  0.0013  -0.0006283  0.0004369 

Table  10.  Comparison  of  Green's  Function  with  LSVOCN  for  Test  Case  3 
(Vary  the  Receiver  Depth)  -  n2(y)  with  Negative  and  Positive  Gradients. 

C.  COMPARISON  OF  LSVOCN  WITH  THE  RECURSIVE  RAY  ACOUSTICS 
ALGORITHM 

1.  Test  Preparation 

The  RRA  Algorithm  calculates  the  magnitude  in  Pascals  and  phase  in  degrees  of 
the  acoustic  pressure  along  a  ray  path.  The  source  is  modeled  as  a  time-harmonic, 
omnidirectional  point  source,  with  its  strength  represented  by  a  source  level  (SL)  value  in 
decibels  relative  to  one  micro-pascal  (rms).  The  "Airy  function  solution"  is  in  terms  of 
velocity  potentials  and  is  based  on  a  time-harmonic,  omnidirectional  point  source 
modeled  by  a  unit  amplitude  impulse  function.  Thus,  before  the  magnitudes  and  phases 
from  the  RRA  Algorithm  and  LSVOCN  (the  "Airy  function  solution")  could  be 
compared,  the  velocity  potential  expressions  were  transformed  into  acoustic  pressure 
expressions  using  the  following  equation  for  time-harmonic  fields  [Ref.  2,  Chapter  2]: 

Pf(r,y)  =  -j2%fp0y/(r,y)  (3.16) 

where/  is  the  frequency  of  the  source  in  Hz,  and  p0  is  the  ambient  or  equilibrium 
density  of  the  fluid  medium  in  kg/m3. 

In  addition,  the  magnitudes  from  LSVOCN  had  to  be  made  equivalent  to  the 
magnitudes  computed  by  the  RRA  Algorithm  by  relating  the  source  level  in  decibels  to 
the  amplitude  of  the  impulse  function  used  to  model  the  omnidirectional  point  source  in 
LSVOCN.  This  was  accomplished  by  deriving  a  magnitude  scale  factor  that  relates  a 
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source  level  in  decibels  to  the  source  strength  of  a  spherical  sound  source  operating  in  the 
monopole  mode  of  vibration  as  the  radius  of  the  sphere  approaches  zero  [see  Eq.(3.3)]. 
Thus,  substituting  Eq.(3.3)  into  Eq.(3.16)  yields 

PfiR , e, V)  =  +j/poSo -eXP^X/?)  ,R>0.  (3.17) 

The  magnitude  of  Eq.(3.17)  is  given  by 

|p/(/?,e,V)|=~5o,il>0.  (3.18) 


Thus,  the  peak  acoustic  pressure,  P0 ,  at  R  equal  to  one  meter  is 

Po  =  M«.e.¥)|fc,  =4£s„. 

Therefore,  the  source  strength,  So ,  in  m3/s,  is  given  by 


(3.19) 


(3.20) 


The  peak  acoustic  pressure,  P0,  in  Pascals  [Ref.  2,  Chapter  1]  due  to  a  given  source  level 
in  decibels  is  given  by: 


(3.21) 


where  P~f  is  the  reference  pressure  equal  to  one  micro-pascal  (rms),  and  SL  is  the  source 
level  of  the  source  in  decibels  relative  to  Pnf.  Since  the  source  level  used  by  the  RRA 
Algorithm  is  a  known  constant  (180.0  dB),  the  peak  acoustic  pressure  can  be  calculated 
using  Eq.(3.21).  This  result  can  then  be  substituted  into  Eq.(3.20)  allowing  the  source 
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strength  to  be  computed.  The  resulting  source  strength  value,  S0,  is  the  magnitude  scale 
factor  that  was  used  to  multiply  the  output  from  LSVOCN. 

In  summary,  the  conversion  process  to  make  the  LSVOCN  magnitude  and  phase 
equivalent  to  that  computed  by  the  RRA  Algorithm  is  as  follows  [multiply  the  right-hand 
side  of  Eq.(3.16)  by  Eq.(3.20)]: 

Pf(r,y)  =  -j2nf  poj^-(pf(r,y)  =  -j4icP0y/(r,y) ,  (3.22) 

where  P0  is  given  by  Eq.(3.21).  Equation  (3.22)  was  added  as  FORTRAN  computer  code 
to  LSVOCN  so  that  the  conversion  process  would  be  taken  care  of  automatically. 

2.  Comparison  of  the  RRA  Algorithm  with  LSVOCN  when  n2(y)  has  a 
Negative  Gradient 

The  RRA  Algorithm  was  used  to  simulate  various  ray  paths  propagated  through 
the  ocean  medium  when  the  square  of  the  index  of  refraction  was  a  linear  function  of 
depth  with  a  negative  gradient,  corresponding  to  a  positive  gradient  sound-speed  profile. 
Two  ray  paths  identified  by  the  initial  launch  angle,  j5c ,  were  selected  for  each  of  the  test 
frequencies,  1000  Hz,  250  Hz,  and  50  Hz.  The  sound-speed  profile  and  the  ray 
propagation  path  for  each  ray  is  illustrated  in  Figure  10.  Launch  angles  were  selected  so 
that  the  ray  paths  did  not  interact  with  the  ocean  surface  or  bottom,  thus  ensuring  that  the 
magnitude  and  phase  calculated  by  the  RRA  Algorithm  were  for  a  free-space  propagation 
problem.  Test  points  along  each  ray  path  were  selected.  The  depth  (receiver  depth)  and 
the  horizontal  range  (HRNG)  corresponding  to  these  points  were  recorded  and  used  as 
inputs  to  LSVOCN.  The  comparison  between  the  RRA  Algorithm  and  LSVOCN 
magnitude  and  phase  values  are  summarized  in  Tables  1 1  through  13,  where  the  asterisk 
means  that  the  entire  row  corresponds  to  data  at  the  location  of  a  turning  point.  The 
calculated  values  in  Tables  1 1  through  13  were  obtained  as  follows: 

1 .  Delta  magnitude  =  RRA  Magnitude  -  LSVOCN  Magnitude 

2.  Delta  phase  =  RRA  Phase  -  LSVOCN  Phase 
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3.  Initial  phase  offset  —  Delta  phase  of  the  first  point  for  each  ray 

4.  Adj  RRA  Phase  =  (RRA  Phase  -  Initial  phase  offset) 

a)  Add  360°  if  Adj  RRA  Phase  angle  is  negative 

b)  Subtract  360°  if  Adj  RRA  Phase  angle  is  greater  than  360° 

c)  Do  not  add  or  subtract  if  Adj  RRA  Phase  angle  is  positive  and  less  than  360° 


5.  %  Difference  Magnitude  = 


RRA  Magnitude  -  LSVOCN  Magnitude 
LSVOCN  Magnitude 


x  100% 


6.  %  Difference  Phase  = 


Adj  RRA  Phase  -  LSVOCN  Phase 

LSVOCN  Phase  Xl00% 


Analysis  of  the  data  if  Tables  11  through  13  indicate  that  the  magnitudes  and 
phases  calculated  by  the  RRA  Algorithm  match  reasonably  well  with  the  magnitudes  and 
phases  calculated  by  LSVOCN.  Note  that  there  is  an  initial  phase  difference  that  exists 
between  the  RRA  Algorithm  and  LSVOCN.  The  RRA  Algorithm  begins  with  zero  phase 
and  calculates  the  phase  thereafter  based  on  travel  time  as  a  ray  propagates.  In  order  to 
correct  for  this  initial  phase  difference,  each  RRA  phase  angle  after  the  first  data  point  on 
each  ray  is  adjusted  by  an  amount  equal  to  the  delta  phase  of  the  first  data  point.  After 
the  RRA  phase  has  been  adjusted,  a  valid  comparison  between  the  RRA  phase  and 
LSVOCN  phase  can  be  performed.  It  can  further  be  observed  that  the  percentage  phase 
error  tends  to  increase  as  the  ray  propagates.  This  seems  reasonable  since  small  errors 
made  in  calculating  travel  time  and,  hence,  phase  in  the  beginning  accumulate  as  the  ray 
propagates  further  from  the  source.  These  phase  errors  could  possibly  be  reduced  by 
selecting  a  smaller  arc  length  step  size  used  by  the  RRA  Algorithm. 
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Figure  10.  Raytrace  when  the  square  of  the  index  of  refraction  is  a  linear  function  of 
depth  with  a  negative  gradient  for  two  separate  ray  paths  with  f}0  =  82°  and  90°. 
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SURFACE  SOUND-SPEED  =  1500  m/s,  SOURCE-SOUND  SPEED  =  1515  m/s 
SOURCE  DEPTH  =  150  m,  n2(y)  WITH  NEGATIVE  GRADIENT 


Receiver 

Receiver  RRA 

RRA  LSVOCN 

LSVOCN 

DELTA 

DELTA 

Adj  RRA 

%  DIFF 

Depth 

HRNG  Mag 

Phase  Mag 

Phase 

Mag 

Phase 

Phase 

Mag 

meters 

meters  Pascals 

Degree  Pascals 

Degrees 

Pascals 

Degrees 

Degrees 

%  DIFF 
Phase 


150.139 


211.757 


256.408 


283.987 


0.99 


500.178 


1000.164 


1500.38 


*  294.58  I  2056.996 


223.48  I  3499.913 


4499.953 


4988.762 


1414.21 


2.80646 


1.40685 


0.94006 


0.68751 


0.40689 


0.31804 


0.28757 


122.377 


17.022 


215.202 


219.763 


218.93 


145.679 


348.299 


238.7271 


1414.06 


2.79558 


1.39683 


0.9673 


0.6809 


0.43452 


0.2707 


0.30947 


212.411 


0.15  -90.034 


107.978  1  0.01088  -90.956 


303.34  0.01002  -88.138 


307.452  -0.02725  -87.689 


304.917  |  0.00661  -85.987 


-0.02762  -95.521 


0.04733  272.413 


337.044  I  -0.0219  |  -98.317 


N/A 


107.056 


305.236 


309.797 


308.964 


235.713 


78.333 


326.865 


0.01061 


0.38919 


0.71734 


-2.81609 


0.97077 


-6.35874 


17.48799 


-7.07661 


N/A 


-0.85388 


0.62504 


0.76272 


1.32725 


-2.27488 


3.22457 


-3.02008 


1414.21  122.376  1410,77  212.154  3.44  -89.778 

J41.636  499.907  2.82896  325.936  2.83278  56.813  ~00382  269.127 

J16.533  999.753  1.41456  39.786  1,44011  131.49  -0,02555  -91.704 

74.652  1499.98  0.94282  209.696  0.96954  302.063  -0.02672  -92.367  299,474 

! 6^029  1 2000.033  I  0.7071  341.127  0.66581  69.119  0.04129  272.008  70.905 


N/A 


55.714 


129.564 


0.24384  I  N/A 


-0.13485  -1.93442 


-1.77417  -1.46475 


-2.75595  -0.85711 


6.20147 


Table  1 1  Comparison  of  RRA  Algorithm  with  LSVOCN  Magnitude  and  Phase  Values  for  different 
Receiver  Depth  and  Horizontal  Range  Points  Located  on  Two  Separate  Sound  Rays,  Identified  by  $0 


SURFACE  SOUND-SPEED  =  1500  m/s,  SOURCE  SOUND-SPEED  =  1515  m/s 
SOURCE  DEPTH  =  150  m,  n2(y)  WITH  NEGATIVE  GRADIENT 
_  FREQUENCY  -  250  HZ,  TURNING  POINT  DATA _ 

Receiver  Receiver  RRA  ^RA  LSVOCN  I  LSVOCN  DELTA  DELTA  Adj  RRA  %DIFF  o/oDIfF 

Depth  HRNG  Mag  Phase  Mag  Phase  Mag  Phase  Phase  Mag  Phase 


Degrees 


Degrees  Degrees 


150.139 


211.757 


256.408 


283.987 


0.99  1414.21 


500.178  2.80646 


1000.164  1.40685 


1500.38  0.94006 


*  294.58  2056.996  0.68751 


223.48  13499.913  0.40689 


4499.953  0.31804 


4988.7621  0.28757 


150  I  1  I  1414.21 


141.636  499.907  I  2.82896 


116.533  I  999.753  1.41456 


74.652  |  1499.98  0.94282 


16.029  12000.0331  0.7071 


300.594 


4.255 


233.8 


234.941 


234.732 


36.42 


267.075 


59.682 


300.594 


171.484 


279.946 


232.424 


355.282 


2.80952 


1.37561 


0.90205 


0.65279 


0.41643 


0.27485 


0.2389 


IU1U21I 


2.82851 


1.42205 


0.98967 


0.7486 


30.646 


93.273 


325.296 


326.776 


321.869 


133.076 


352.622 


151.482 


-0.86 


-0.00306 


0.03124 


0.03801 


0.03472 


0.04319 


0.04867 


30.689 


260.489  I  0.00045 


11.855  -0.00749 


322.027  -0.04685 


87.039  !  -0.0415 


269.948 


-89.018 


-91.496 


-91.835 


-87.137 


-96.656 


-85.547 


-91.8 


269.905 


-89.005 


268.091 


-89.603 


268.243 


N/A 


94.307 


323.852 


324.993 


324.784 


126.472 


357.127 


149.734 


-0.06077 


-0.10892 


2.27099 


4.21374 


5.31871 


-2.2909 


15.71403 


N/A 


1.10857 


-0.4439 


-0.54563 


0.90565 


-4.96258 


1.27757 


20.37254  1  -1.15393 


N/A  J  -0.09396  I  N/A 


261.579  |  0.01591  I  0.41844 


10.041  |  -0.5267  -15.30156 


322.519  |  -4.7339  0.15278 


85.377  |  -5.54368  I  -1.90949 


Table  12.  Comparison  of  RRA  Algorithm  with  LSVOCN  Magnitude  and  Phase  Values  for  different 
Receiver  Depth  and  Horizontal  Range  Points  Located  on  Two  Separate  Sound  Rays,  Identified  by  p0 
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SURFACE  SOUND-SPEED  =  1500  m/s,  SOURCE  SOUND-SPEED  =  1515  m/s 
SOURCE  DEPTH  =  150  m,  n2(y)  WITH  NEGATIVE  GRADIENT 
FREQUENCY  =  50  HZ,  TURNING  POINT  DATA 


Receiver  Receiver  RRA  RRA  LSVOCN  LSVOCN  DELTA 

Depth  HRNG  Mag  Phase  Mag  Phase  Mag 

meters  meters  Pascals  Degree  Pascals  Degrees  Pascals 


%DIFF 

%DIFF 

Mag 

Phase 

150.139 


211.757 


256.408 


283.987 


♦  294.58 


223.48 


90.742 


1 


150 


141.636 


116.533 


74.652 


16.029 


0.99 


500.178 


1000.164 


1500.38 


2056.996 


3499.913 


4499.953 


4988.762 


1414.21 


2.80646 


1.40685 


0.94006 


0.68751 


0.40689 


0.31804 


0.28757 


348.119 


144.851 


334.76 


190.988 


118.946 


7.284 


269.415 


83.936 


1396.42 


2.7779 


1.44067 


0.89268 


0.7305 


0..3614 66 


0.36089 


0.32794 


77.964 


235.583 


63.371 


281.376 


207.172 


99.872 


17.79 


0.02856 


-0.03382 


0.04738 


-0.04299 


0.04542 


168.751  -0.04037 


270.155 


-90.732 


271.389 


-90.388 


-88.226 


-92.588 


-86.008 


-84.815 


N/A 


234.696 


64.605 


280.833 


208.791 


97.129 


359.26 


173.781 


1.27397 


1.02811 


-2.34752 


5.30761 


-5.88501 


12.66604 


-11.87342 


-12.31018 


0.78148 


-2.74652 


1.07956 


2.98072 


499.907 


999.753 


1499.98 


2000.033 


1414.21 


2.82896 


1.41456 


0.94282 


0.7071 


348.119 


178.297 


343.989 


118.485 


215.056 


1425.34 


2.8749 


1.4166 


0.95694 


0.67256 


78.211 


267.929 


75.872 


205.742 


302.346 


^^ME£Ei31 
EmiiMIEI 


-11.13 


0.00204 


-0.01412 


0.03454 


N/A 

-0.78087 

268.389 

-1.59797 

74.081 

-0.14401 

208.577 

-1.47554 

305.148 

5.1356 

0.17169 


-2.36055 


1.37794 


0.92675 


Table  13.  Comparison  of  RRA  Algorithm  with  LSVOCN  Magnitude  and  Phase  Values  for  different 
Receiver  Depth  and  Horizontal  Range  Points  Located  on  Two  Separate  Sound  Rays,  Identified  by  $0 


3.  Comparison  of  the  RRA  Algorithm  with  LSVOCN  when  n2(y)  has  a 
Positive  Gradient 

The  RRA  Algorithm  was  used  to  simulate  various  ray  paths  propagated  through 
the  ocean  medium  when  the  square  of  the  index  of  refraction  was  a  linear  function  of 
depth  with  a  positive  gradient,  corresponding  to  a  negative  gradient  sound-speed  profile. 
Two  ray  paths  identified  by  the  initial  launch  angle,  (JQ ,  were  selected  for  each  of  the  test 
frequencies,  1000  Hz,  250  Hz,  and  50  Hz.  The  sound-speed  profile  and  the  ray 
propagation  path  for  each  ray  is  illustrated  in  Figure  1 1 .  Launch  angles  were  selected  so 
that  the  ray  paths  did  not  interact  with  the  ocean  surface  or  bottom,  thus  ensuring  that  the 
magnitude  and  phase  calculated  by  the  RRA  Algorithm  were  for  a  free-space  propagation 
problem.  Test  points  along  each  ray  path  were  selected.  The  depth  (receiver  depth)  and 
the  horizontal  range  (HRNG)  corresponding  to  these  points  were  recorded  and  used  as 
inputs  to  LSVOCN.  The  comparison  between  the  RRA  Algorithm  and  LSVOCN 
magnitude  and  phase  values  are  summarized  in  Tables  14  through  16,  where  the  asterisk 
means  that  the  entire  row  corresponds  to  data  at  the  location  of  a  turning  point.  The 
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calculated  values  in  Tables  14  through  16  were  obtained  in  the  same  manner  as  was  done 
for  Tables  1 1  through  13.  Analysis  of  the  data  in  Tables  14  through  16  indicate  that  the 
magnitudes  and  phases  calculated  by  the  RRA  Algorithm  match  reasonably  well  with  the 
magnitudes  and  phases  calculated  by  LSVOCN.  Note  that  there  is  an  initial  phase 
difference  that  exists  between  the  RRA  Algorithm  and  LSVOCN.  The  RRA  Algorithm 
begins  with  zero  phase  and  calculates  the  phase  thereafter  based  on  travel  time  as  a  ray 
propagates.  In  order  to  correct  for  this  initial  phase  difference,  each  RRA  phase  angle 
after  the  first  data  point  on  each  ray  is  adjusted  by  an  amount  equal  to  the  first  delta 
phase  of  the  first  data  point.  After  the  RRA  phase  has  been  adjusted,  a  valid  comparison 
between  the  RRA  phase  and  LSVOCN  phase  can  be  performed.  It  can  further  be 
observed  that  the  percentage  phase  error  tends  to  increase  as  the  ray  propagates.  This 
seems  reasonable  since  small  errors  in  the  beginning  accumulate  as  the  ray  propagates 
further  from  the  source.  These  phase  errors  could  possibly  be  reduced  by  selecting  a 
smaller  arc  length  step  size  used  by  the  RRA  Algorithm.  Also,  for  small  acoustic 
pressure  phases,  small  errors  result  in  large  percentage  differences,  as  expected. 
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Table  14.  Comparison  of  RRA  Algorithm  with  LSVOCN  Magnitude  and  Phase  Values  for  different 
Receiver  Depth  and  Horizontal  Range  Points  Located  on  Two  Separate  Sound  Rays,  Identified  by  (}0 . 
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SURFACE  SOUND-SPEED  =  1500  m/s,  SOURCE  SOUND-SPEED  =  1485 
SOURCE  DEPTH  =  150  m,  n2(y)  WITH  POSITIVE  GRADIENT 
FREQUENCY  =  250  HZ,  TURNING  POINT  DATA 


LSVOCN 

Mag 

Pascals 


Receiver 

Receiver 

RRA 

Depth 

HRNG 

Mag 

meters 

meters 

Pascals 

DELTA 

DELTA 

Mag 

Phase 

Pascals 

Decrees 

%DIFF 

%DIFF 

Mag 

Phase 

149.861 


88.159 


43.255 


15.25 


*  3.967 


12251 


202.264 


290.462 


p=90° 


158.268 


183.11 


224.566 


282.604 


0.99 


500.168 


1000.132 


1500.325 


2077.678 


3498.868 


4499.27 


4988.36 


499.534 


999.394 


1499.657 


1999.778 


1414.21 


2.80645 


1.40682 


0.94002 


0.68067 


0.40695 


0.31802 


0.28753 


1414.21 


2.83107 


1.41507 


0.94303 


0.70719 


299.394 


120.102 


109.303 


232.708 


124.586 


247.665 


10.566 


300.064 


299.394 


314.164 


182.002 


252.073 


130.724 


1415.07 


2.77775 


1.44734 


0.98028 


0.67268 


0.4037 


0.30859 


0.33218 


1406.07 


2.80951 


1.37475 


0.94032 


0.67483 


29.447 


209.334 


200.284 


321.198 


218.573 


344.465 


91.865 


33.548 


29.774 


44.984 


272.98 


339.194 


223.51 


269.947 


0.0287  -89.232 


-0.04052  -90.981 


-0.04026  -88.49 


0.00799 


0.00325 


0.00943 


-0.04464 


8.14 


0.02156 


0.04032 


0.00271 


0.03236 


-93.987 


-96.8 


-81.299 


266.516 


269.62 


269.18 


-90.978 


-87.121 


-92.786 


N/A 


210.155 


199.356 


322.761 


214.639 


337.718 


100.619 


30.117 


N/A 


44.544 


272.382 


342.453 


221.104 


-0.06077 


1.03321 


-2.79962 


—4.10699 


1.18779 


0.80505 


3.05583 


-13.44151 


0.57892 


0.76739 


2.9329 


0.2882 


4.79528 


N/A 


0.3922 


-0.46334 


0.48662 


-1.79986 


-1.95869 


9.5292 


N/A 


-0.97813 


-0.21906 


0.96081 


-1.07646 


Table  15.  Comparison  of  RRA  Algorithm  with  LSVOCN  Magnitude  and  Phase  Values  for  different 
Receiver  Depth  and  Horizontal  Range  Points  Located  on  Two  Separate  Sound  Rays,  Identified  by  . 


SURFACE  SOUND-SPEED  -  1500  m/s,  SOURCE  SOUND-SPEED  =  1485  m/s 
SOURCE  DEPTH  =  150  m,  n2(y)  WITH  POSITIVE  GRADIENT 
FREQUENCY  =  50  HZ,  TURNING  POINT  DATA 


DELTA 

Mag 

Pascals 


Receiver 

Receiver 

Depth 

HRNG 

meters 

meters 

%DIFF 

%DIFF 

Mag 

Phase 

149.861 


88.159 


43.255 


15.25 


*  3.967 


72.257 


202.264 


290.462 


Po=90° 


150 


158.268 


183.11 


224.566 


282.604 


0.99 


500.168 


1000.132 


1500.325 


2077.678 


3498.868 


4499.27 


4988.36 


499.534 


999.394 


1499.657 


1999.778 


1414.21 


2.80645 


1.40682 


0.94002 


0.68067 


0.40695 


0.31802 


0.28753 


1414.21 


2.83107 


1.41507 


0.94303 


0.70719 


347.879 


24.02 


93.861 


190.542 


96.917 


265.533 


290.113 


348.013 


347.879 


62.833 


108.4 


122.415 


98.145 


1395.57 


2.8002 


1.37231 


0.89886 


0.72557 


0.43785 


0.36357 


0.3022 


1417.16 


2.82662 


1.45024 


0.94316 


0.74101 


\msm\ 


182.489 


279.047 


188.197 


350.513 


17.357 


86.994 


78.123 


153.724 


197.13 


209.544 


185.56 


18.64 


0.00625 


0.03451 


0.04116 


-0.0449 


-0.04555 


-0.01467 


-2.95 


0.00445 


-0.03517 


-0.00013 


270.166 


-88.966 


272.756 


261.019 


269.756 


-90.891 


msm* 


183.695 


280.376 


186.751 


355.367 


19.947 


77.847 


1.33565 


0.2232 


2.51474 


4.57913 


-6.18824 


-7.05721 


-12.52854 


-4.8544 


N/A 

-0.20816 

153.077 

0.15743 

198.644 

-2.42512 

212.659 

-0.01378 

188.389 

-4.56404 

0.66086 


0.47626 


-0.76834 


1.38483 


14.92193 


-10.51452 


0.76802 


1.48656 


1.52457 


Table  16.  Comparison  of  RRA  Algorithm  with  LSVOCN  Magnitude  and  Phase  Values  for  different 
Receiver  Depth  and  Horizontal  Range  Points  Located  on  Two  Separate  Sound  Rays,  Identified  by  . 
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Figure  11.  Raytrace  when  the  square  of  the  index  of  refraction  is  a  linear  function  of 
depth  with  a  positive  gradient  for  two  separate  ray  paths  with  po  =  90°  and  (}„  =  98°. 
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IV.  SUMMARY  AND  RECOMMENDATIONS 

This  thesis  had  two  primary  goals.  First,  derive  and  document  the  solution  to  the 
three-dimensional  inhomogeneous  Helmholtz  equation  for  a  ffee-space  propagation 
problem  when  the  square  of  the  index  of  refraction  of  the  ocean  medium  is  a  linear 
function  of  depth,  and  the  source  is  an  omnidirectional  point  source.  Given  this  type  of 
index  of  refraction,  the  solution  to  the  inhomogeneous  Helmholtz  equation  has  an  exact 
solution  in  terms  of  Airy  functions.  This  exact  solution,  in  terms  of  Airy  functions,  was 
then  incorporated  into  the  computer  program  Linear  Space-Variant  Ocean  (LSVOCN)  as 
an  additional  ocean  medium  transfer  function.  The  second  goal  was  to  further  test  the 
Recursive  Ray  Acoustics  (RRA)  Algorithm  by  comparing  the  magnitudes  and  phases  of 
the  acoustic  sound  pressure  calculated,  at  various  points,  along  a  ray  path  generated  by 
the  RRA  Algorithm  with  corresponding  magnitudes  and  phases  calculated  by  LSVOCN 
along  the  same  ray  path. 

In  Chapter  II,  the  solution  to  the  aforementioned  Helmholtz  equation  was 
carefully  derived  and  documented.  This  documentation  served  to  record  in  detail  the 
steps  that  led  to  the  final  Airy  function  solution. 

The  next  step  in  the  solution  process  was  to  convert  the  theoretical  Airy  function 
solution  to  the  FORTRAN  computer  code  that  was  used  by  LSVOCN.  The  main 
problems  in  the  conversion  process  were  the  many  numerical  errors  that  developed. 

There  were  two  main  numerical  errors  encountered,  overflow  and  underflow  errors. 
Overflow  errors  resulted  from  the  evaluation  of  the  Bairy  function  with  large  positive 
argument  values.  The  Airy  function  did  not  have  this  problem  since  it  is  well  behaved  for 
both  negative  and  positive  argument  values.  A  graphical  representation  of  both  the  Airy 
and  Bairy  functions  is  shown  in  Figures  7  and  8,  respectively.  Overflow  errors  result  in 
the  termination  of  the  computer  code  before  it  has  finished  a  calculation.  Thus  no  answer 
is  ever  obtained  -  a  serious  problem.  To  prevent  overflow  errors,  we  used  the 
exponentially  scaled  Airy  and  Bairy  functions,  also  shown  in  Figures  7  and  8.  These  two 
functions  allowed  a  significant  advance  towards  a  working  FORTRAN  computer  code 
that  would  solve  the  "Airy  function  solution"  as  given  by  Eq.(2.75)  through  Eq.(2.78). 


49 


Using  the  exponentially  scaled  functions  allowed  for  very  large  positive  values  of  £(y), 
the  argument  of  the  Airy  and  Bairy  functions,  to  be  used.  These  large  positive  values  of 
£(y)  were  the  result  of  propagating  waves  becoming  evanescent  waves.  For  the  solution 
to  be  accurate,  it  is  very  important  to  account  for  these  evanescent  waves.  The 
exponentially  scaled  functions  did  not  prevent  underflow  errors.  These  types  of  errors  do 
introduce  some  numerical  inaccuracies,  but  do  not  result  in  computer  program 
termination. 

The  last  numerical  problem  to  solve  was  to  establish  an  upper  limit  of  integration 
for  Eq.(2.3).  Using  the  computer,  it  is  impossible  to  literally  integrate  to  infinity.  To 
achieve  numerical  accuracy,  the  upper  limit  of  integration  had  to  be  large  enough  to 
account  for  most  of  the  significant  evanescent  waves.  Thus,  a  method  to  automatically 
compute  this  upper  limit  of  integration,  regardless  of  the  source  frequency  selected,  was 
developed. 

Before  comparing  the  RRA  Algorithm  with  LSVOCN,  we  developed  test  cases  to 
validate  the  Airy  function  solution  and  FORTRAN  computer  code  used  by  LSVOCN. 
This  validation  was  carried  out  by  comparing  the  magnitude  and  phase  of  the  velocity 
potential  equal  to  the  free-space  Green’s  function  of  a  homogeneous  ocean  medium  with 
the  magnitude  and  phase  of  the  velocity  potential  computed  by  LSVOCN.  The 
free-space  Green's  function  was  chosen  because  we  could  make  the  modeled  ocean 
medium  used  by  LSVOCN  "simulate"  a  homogeneous  medium  without  actually  being 
made  homogeneous.  Therefore,  the  Green’s  function  provided  a  simple  theoretical 
solution  that  we  could  compare  the  calculated  solution  of  LSVOCN  with.  All  three 
Green's  function  test  cases  discussed  in  detail  in  Chapter  III  gave  very  good  results.  Now 
we  were  finally  ready  to  test  calculated  values  from  the  Recursive  Ray  Acoustics  (RRA) 
Algorithm  with  those  predicted  by  LSVOCN. 

The  RRA  Algorithm  versus  LSVOCN  comparison  tests  were  conducted  first, 
when  the  square  of  the  index  of  refraction  had  a  negative  gradient  and  second,  when  the 
square  of  the  index  of  refraction  had  a  positive  gradient.  The  test  results  indicated  that 
the  magnitudes  from  the  RRA  Algorithm  and  LSVOCN  were  in  close  agreement.  In 
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order  to  be  able  to  compare  the  RRA  Algorithm  with  LSVOCN,  a  magnitude  scale  factor 
was  derived  that  related  the  source  level  in  decibels  used  by  the  RRA  Algorithm  to  the 
amplitude  of  the  impulse  function  used  to  model  the  omnidirectional  point  source  in 
LSVOCN.  The  magnitudes  of  LSVOCN  were  extremely  close  to  the  values  predicted  by 
the  RRA  Algorithm. 

The  phase  comparisons  between  the  RRA  Algorithm  and  LSVOCN  were  also 
quite  good.  The  initial  phase  offset  between  LSVOCN  and  the  RRA  Algorithm  was 
determined  by  calculating  the  phase  differences,  one  meter  from  the  source,  between  the 
RRA  Algorithm  and  LSVOCN,  i.e.,  the  phase  difference  between  the  RRA  Algorithm 
and  LSVOCN  of  the  first  data  point  of  each  ray  path.  Once  this  difference  was 
determined,  it  was  subtracted  from  the  other  phase  values,  less  the  first  data  point, 
calculated  by  the  RRA  Algorithm  along  the  same  ray  path.  It  was  observed  that  the 
percentage  phase  error  tended  to  grow  as  the  receiver  was  moved  further  from  the  source. 
This  trend  makes  intuitive  sense  because  it  suggests  that  small  phase  errors  in  the 
beginning  are  propagated  and  accumulated  as  the  ray  travels  further  from  the  source.  It 
was  also  observed  that  several  data  points  resulted  in  large  percentage  differences  for  the 
phase.  These  large  percentage  differences,  resulting  from  the  variance  between  small 
phases  could  possibly  be  reduced  by  making  the  integration  step  size  a  function  of  the 
wave  number  rather  than  a  fixed  quantity. 

Because  ray  acoustics  is  an  approximate  solution  of  the  wave  equation, 
characterized  as  a  "high  frequency"  approximation,  we  expected  to  see  better  and  better 
agreement  between  the  RRA  Algorithm  and  LSVOCN  as  the  source  frequency  increased 
from  50  Hertz  to  1000  Hertz.  The  data  in  Tables  1 1  through  1 6  did  not  indicate  such  a 
trend  for  all  cases.  In  fact,  some  percentage  errors  in  magnitude  and  phase  are  larger  for 
the  1000  Hertz  case  than  the  50  Hertz  case.  The  expected  results  could  possibly  be 
observed  if  the  integration  step  size  were  made  a  function  of  the  wave  number,  as 
discussed  earlier.  The  results  from  all  six  test  cases  comparing  the  RRA  Algorithm  with 
LSVOCN  indicate  that  the  RRA  Algorithm  was,  in  general,  accurate  regardless  of  the 
frequency  selected. 
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As  a  final  note,  the  CPU  time  required  by  the  RRA  Algorithm  to  compute  all  of 
the  magnitudes  and  phases  for  thousands  of  points  along  each  ray  path  was  typically  on 
the  order  of  30  seconds.  Table  17  gives  the  approximate  CPU  times  for  each  LSVOCN 
data  point  calculated.  Examination  of  Table  17  reveals  that  the  RRA  Algorithm  is  a  much 
more  efficient  method  to  calculate  the  acoustic  pressure  compared  to  the  Airy  function 
wave  solution.  The  wave  solution  not  only  takes  considerably  more  time,  but  is  also  more 
susceptible  to  numerical  inaccuracies. 


CPU  TIMES  PER  DATA  POINT  CALCULATED  BY  LSVOCN  | 

Po=  90° 

1  P<r  82° 

P<r  90° 

[  P0=  98° 

Square  of  the  index  of  refraction  with 
Negative  gradient 

Square  of  the  index  of  refraction  with 
Positive  gradient 

receiver  @  1  meter 
~59  min 

other  receiver  points 
~6  min 

receiver  @  1  meter 
~50  min 

other  receiver  points 
~6  min 

receiver  @  1  meter 
~54  min 

other  receiver  points 
~6  min 

receiver  @  1  meter 
~22.5  min 

other  receiver  points 
~6  min 

■ 

receiver  @  1  meter 
~  16  min 

other  receiver  point 
~5  min 

receiver  @  1  meter 
~  22  min 

other  receiver  point 
~5  min 

receiver  @  1  meter 
~  16  min 

other  receiver  point 
~5.5  min 

receiver  @  1  meter 
~  23  min 

other  receiver  point 
~5.3  min 

■ 

receiver  @  1  meter 
~8  min 

other  receiver  point 
~5  min 

receiver  @  1  meter 
~25  min 

other  receiver  point 
~5  min 

receiver  @  1  meter 
~8.2  min 

other  receiver  point 
~4  min 

receiver  @  1  meter 
~26  min 

other  receiver  point 
~5  min 

Table  17.  CPU  Times  Per  Data  Point  Calculated  by  LSVOCN. 


Recommendations  for  further  research  include: 

1.  To  achieve  more  accurate  calculated  acoustic  pressure  values  from  the  RRA 
Algorithm  (magnitude  and  phase)  a  smaller  arc  length  step  size,  used  by  the 
RRA  Algorithm,  should  be  selected.  Using  a  smaller  arc  length  step  size  may 
decrease  small  phase  errors  that  are  propagated  and  accumulated  as  the  ray 
travels  further  from  the  source. 

2.  Propagate  sound  rays,  using  the  RRA  Algorithm,  such  that  focal  points  are 
generated.  Compare  the  RRA  Algorithm  and  LSVOCN  results. 

3.  Make  the  integration  step  size  used  by  Eq.(2.3)  a  function  of  the  wave  number 
to  determine  if  improved  accuracy  can  be  obtained. 
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APPENDIX.  TABULATED  VALUES  FOR  n2(y),kj(y),an&  £(y) 
FOR  DIFFERENT  LAUNCH  ANGLES  AND  RECEIVER  DEPTHS 


FREQUENCY  =  50  HZ 

SOURCE  DEPTH  =  150  M 
n2(y)  WITH  POSITIVE  GRADIENT 

RECEIVER  DEPTH 

1— 

>  =  50 
meters 

>  =  100 
meters 

y  =  150 
meters 

>  =  200 
meters 

>  =  250 
meters 

>  =  300 
meters 

S  ^  -r* 

COL 

0.98010 

-0.0004016 

1.2249 

0.98673 
-0.0001048 
0.31946  | 

0.99337 
0.00019213 
-0.58594  | 

1.00000 

0.00048901 

-1.4913 

1.0066 

0.00078589 

-2.3967 

1.0133 

0.0010828 

-3.3021 

1.0199 

0.0013796 

-4.075 

P0  =  86° 
■2(y) 
k,J(y) 

C(y> 

0 . 98010 
-0.0006729 
2 . 0520 

0.98673 

-0.0003759 

1.1466 

0.99337 

-0.0000791 

0.24123 

1.00000 

0.00021778 

-0.66417 

1.0066 

0.00051466 

-1.5696 

1.0133 

0.00081154 

-2.4750 

1.0199 

0.0011084 

-3.3804 

P„  =  88° 

■to) 

*,to> 

C(y) 

0.98010 

-0.0008361 

2.5500 

0.98673 

-0.0005393 

1.6446 

0.99337 

-0.0002424 

0.73915 

1 . 00000 

0.00005451 

-0.16624 

1.0066 

0.00035139 

-1.0716 

1.0133 

0.00064827 

-1.9770 

1.0199 

0.00094515 

-2.8824 

Po  =  90° 

■to) 

«S2(y) 

«y> 

0.98010 

-0.0008906 

2.7162 

0.98673 

-0.0005938 

1.8108 

0.99337 

-0.0002969 

0.90540 

1 . 00000 

0.0 

0.0 

1.0066 

0.00029688 

-0.90540 

1.0133 

0.00059376 

-1.8108 

1.0199 

0.00089064 

-2.7162 

P0  —  92° 
■to) 

V(y) 

C(y) 

0 . 98010 
-0.0008361 

2 . 5500 

0.98673 
-0 . 0005393 
1.6446 

0.99337 

-0.0002424 

0.73915 

1.00000 

0.00005451 

-0.16624 

1.0066 

0.00035139 

-1.0716 

1.0133 

0.00064827 

-1.9770 

1.0199 

0.00094515 

-2.8824 

Po  =  94° 

■to) 

k,J(y) 

C(y) 

0.98010 

-0.0006729 

2 . 0520 

0.98673 
-0 . 0003760 
1.1466 

0.99337 

0.00019213 

-0.58594 

1.00000 

0.00048901 

-1.4913 

1.0066 

0.00078589 

-2.3967 

1.0133 

0.0010828 

-3.3021 

1 . 0199 

0.0011084 

-3.3804 

Pc  =  96° 
■to) 
k,2(y) 

C(y) 

0 . 98010 
-0.0004016 
1.2249  * 

0.98673 

-0.0001048 

0.31946 

0.99337 

-0.0001028 

0.031345 

1.00000 

0.00048901 

-1.4913 

1.0066 

0.0010847 

-3.3081 

1.0133 

0.0016844 

-5.1370 

1.0199 

0.0013796 

-4.2075 

P0  =  98° 
■to) 

V(y) 

«y> 

0.98010 

-0.0000238 

0.072457 

0.98673 

0.00027312 

-0.83294 

0.99337 
0.00057000 
-1 . 7383 

1.00000 

0.00086688 

-2.6437 

1.0066 

0.0011638 

-3.5491 

1.0133 

0.0014606 

-4.4545 

1.0199 

0.0017575 

-5.3599 

Table  18.  Values  of  n2(y),  kj(y ),  £(y)  for 
different  launch  angles  p0and  receiver  depths  y. 
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FREQUENCY  =  50  HZ 
SOURCE  DEPTH  =  150  M 
n2(y)  WITH  NEGATIVE  GRADIENT 


RECEIVER  DEPTH 


y  =  100 

meters 


y  =  15© 
meters 


y  =  200 
meters 


y  =  250 
meters 


y  =  300 
meters 


Pc  =  84° 
■*(y) 
V(y) 
C(y> 


1.0201  1.0134  1.0067  1.00000  0.99330  0.98660  0.97990 

0.0013341  0.0010460  0.00075794  0.00046983  0.00018173  -0.0001064  -0.0003945 

-4.1510  -3.2546  -2.3582  -1.4618  -0.56542  0.33097  1.2274 


1.00000  0.99330 

0.00020924  -0.0000789 
-0.65101  0.24537 


P0  =  88 

■to) 

V(y) 

C(y> 


P„  =  90' 
■to) 
•sto) 
C(y) 


Po  =  92‘ 

■to) 

>SJ(y) 

«y) 


P0  =  94 
■to) 
Kto) 

«y) 


Pc  =  96° 

■to) 

V(y) 

C(y) 


Po  -  98° 

■to) 

V(y) 

CCv) 


1.0067 

8  0.00034048 
-1.0593 


1.0201  1.0134  1.0067 

0.00086431  0.00057621  0.00028810 
-2.6892  -1.7928  -0.89639 


0.97990 

-0.0006551 

2.0382 


0.98660  0.97990 

-0.0005238  -0.0008119 
1.6298  2.5262 


0.99330 

-0.0002357 

0.73344 


0.97990 
62  -0.0008643 
2.6892 


0.97990 
38  -0.0008119 
2.5262 


0.99330 

0.00018173 

-0.66542 


0.98660  0.97990 

-0.0003670  -0.0006551 
1.1418  2.0382 


0.97990 

-0.0003945 

1.2274 


1.0067  1.00000 

0.0011210  0.00083289 
-3.4878  -2.5914 


0.98660  0.97990 

8  0.00025668  -0.0000314 
-0.79861  0.097780 


Table  21.  Values  of  n2(y),  k2(y),  £(y)  for 
different  launch  angles  Po  and  receiver  depths 


FREQUENCY  =  250  HZ 
SOURCE  DEPTH  =  150  M 
n2(y)  WITH  NEGATIVE  GRADIENT 


r~ 

RECEIVER  DEPTH 

■csss 

■ 

y  =  loo 
meters 

y  =  150 
meters 

y  =  200 
meters 

y  =  250 
meters 

y  =  300 
meters 

P0  =  84° 
■*(y) 
V(y) 

C(y> 

1.0201 

0 . 033654 
-12.138 

1.0134 

0.026151 

-9.5165 

1.0067 

0.018948 

-6.8954 

1 . 00000 

0.011746 

-4.2744 

0.99330 
0.0045432 
-1 . 6533 

0.98660 

-0.0026594 

0.96776 

0.97990 

-0.0098620 

3.5888 

P0  =  86° 
■*00 

«y> 

1.0201 

0.026839 

-9.7667 

1.0134 

0.019636 

-7.1457 

1.0067 

0.012434 

-4.5246 

1.00000 

0.0052310 

-1.9036 

0.99330 

-0.0019716 

0.71748 

0.98660 

-0.0091742 

3.3385 

0.97990 

-0.016377 

5.9596 

Pc  =  88° 

■%) 

V(y) 

C(y> 

1.0201 

0.022917 

-8.3396 

1.0134 

0.015715 

-5.7186 

1.0067 

0.0085119 

-3.0975 

1.00000 

0.0013093 

-0.47647 

0 . 99330 
-0 . 0058933 
2.1446 

0.98660 

-0.013096 

4 . 7656 

0.97990 

-0.020298 

7.3867 

Po  =  90° 
■*(y) 

V(y) 

C(y) 

1.0201 

0.021608 

-7.8632 

1.0134 

0.014405 

-5.2421 

1.0067 

0.0072026 

-2.6211 

0.99330 

-0.0072026 

2.6211 

0.98660 

-0.014405 

5.2421 

0.97990 

-0.021608 

7.8632 

P0  =  92° 
■*(y) 

V(y) 

?(y> 

1.0201 

0.022917 

-8.3396 

1.0134 

0.015715 

-5.7186 

1.0067 

0.0085119 

-3.0975 

1 . 00000 

0.0013093 

-0.47647 

0 . 99330 

-0.0058933 

2.1446 

0.98660 

-0.013096 

4.7656 

0.97990 

-0.020298 

7.3867 

P0  =  94° 
■*(y) 
k,J(y) 

C(y> 

1.0201 

0.026839 

-9.7667 

1.0134 

0.019636 

-7.1457 

1.0067 

0.012434 

-4.5246 

1.00000 
0.0052310 
-1 . 9036 

0 . 99330 

-0.0019716 

0.71748 

0.98660 

-0.0091742 

3 . 3385 

0.97990 

-0.016377 

5.9596 

P0  =  96° 
■ty) 

V(y) 

C(y) 

1.0201 

0 . 033354 
-12.138 

1.0134 

0.026151 

-9.5165 

1.0067 

0.018948 

-6.8954 

1 . 00000 

0.011746 

-4.2744 

0.99330 

0.0045432 
-1 . 6533 

0.98660 

-0.0026594 

0.96776 

0.97990 

-0.0098620 

3.5888 

p0  =  98° 

"’(y) 

t,2(y) 

C(y> 

1.0201 

0.042430 

-15.440 

1.0134 

0.035227 

-12.819 

1.0067 

0 . 028025 
-10.198 

1.00000 

0.020822 

-7.5773 

0.99330 

0.013620 

-4.9562 

0.98660 

0.0064169 

-2.3351 

0.97990 

-0.0007857 

0.28591 

Table  22.  Values  of  n2(y),  k2(y),  £(y)  for 
different  launch  angles  (Jo  and  receiver  depths^. 
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FREQUENCY  =  1000  HZ 
SOURCE  DEPTH  =  150  M 
n2(y)  WITH  NEGATIVE  GRADIENT 


RECEIVER  DEPTH 


y  =  200 
meters 


y  =  250 
meters 


y  =  300 
meters 


P0  =  88° 

■ty) 

V<y) 

C(y) 


P„  =  90° 

^(y) 

*S2(y) 

C(y) 


p0  =  92° 
•^(y) 
«S2(y) 
C(y> 


Po  =  94° 
■ty) 
*S2(y) 
C(» 


Po  =  96° 

n2(y) 

k,2(y) 

CCv) 


Pc  =  98° 

k,2(y) 

C(y> 


0.98660  0.97990 

*0.042550  -0.15779 

2.4386  9.0432 


0.99330  0.98660  0.97990 

-0.031546  -0.14679  -0.26203 

1.8079  8.4126  15.017 


0.97990 

-0.32478 

18.613 


0.97990 

-0.34572 

19.614 


0.97990 

-0.32478 

18.613 


0.97990 

-0.26203 

15.017 


0.98660  0.97990 

-0.042550  -0.15779 

2.4386  9.0432 


1.00000 

0.33315 

-19.094 


0.99330 

0.21791 

-12.489 


0.98660 

0.10267 

-5.8842 


0.97990 

-0.012571 

0.72045 


Table  23.  Values  of  n2(y),  kj(y ),  £(y)  for 
different  launch  angles  Po  and  receiver  depths  j'. 


DEPTH  DEPTH 


0  8  "5  0  6  -6  0  6  -10  0  10 

ZETA  ZETA  ZETA  ZETA 


FREQUENCY  =60  hertz 

SOURCE  DEPTH  =  160  meters 

OCEAN  DEPTH  =  300  meters 

SOUND  SPEED,  OCEAN  SURFACE  =  1800  mfs 

SOUND  SPEED,  SOURCE  DEPTH  =  1486  nYs 


Figure  12.  Plot  of  C(y)  versus  depth  y  for  different  launch  angles  (30 : 
nz(y)  with  positive  gradient. 


ZETA  ZETA  ZETA  ZETA 

FREQUENCY  =  260  hertz 

SOURCE  DEPTH  =  150  meters 

OCEAN  DEFTH  =  300  meters 

SOUND  SPEED,  OCEAN  SURFACE  =  1600  rrYs 

SOUND  SPEED,  SOURCE  DEPTH  =  1486  rrVs 


Figure  13.  Plot  of  £(y)  versus  depth  y  for  different  launch  angles  p0 : 
nz(y)  with  positive  gradient . 
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ZETA  ZETA  ZETA  ZETA 


FREQUENCY  =  1000  hertz 

SOURCE  DEFTH  =  150  meters 

OCEAN  DEPTH  =  900  meters 

SOUND  SPEED,  OCEAN  SURFACE  =  1500  rrVs 

_ _ SOUND  SPEED,  SOURCE  DEPTH  =  1486  rrVs 

Figure  14.  Plot  of  £(y)  versus  depth  y  for  different  launch  angles  (Jo 
n2(y)  with  positive  gradient. 


ZETA  ZETA  ZETA  ZETA 


FREQUENCY  =  50  hertz 
SOURCE  DEPTH  =  150  meters 
OCEAN  DEPTH  =  900  meters 
SOUND  SPEED.  OCEAN  SURFACE  =  1500  nVs 
_ SOUND  SPEED.  SOURCE  DEFTH  =  1616  nVs 

Figure  15.  Plot  of  £(y)  versus  depth  y  for  different  launch  angles  f$0  : 
n2(y)  with  negative  gradient. 
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DEPTH  DEPTH  I  ^  |  DEPTH  DEPTH 


Figure  16.  Plot  of  £(y)  versus  depth  j?  for  different  launch  angles  JJo : 
n2(y)  with  negative  gradient. 


BET AO  =  84  BET AO  =  86  BETAO  =  B8  BETAO  =  90 


-60  0  60  -60  0  60  -50  0  50  -20  0  20 

ZETA  ZETA  ZETA  ZETA 

BETAO  =  92  BETAO  =  94  BETAO  =  96  BETAO  =90 


ZETA  ZETA  ZETA  ZETA 

FREQUENCY  =  1000  hertz 

SOURCE  DEPTH  =  150  meters 

OCEAN  DEPTH  =  300  meters 

SOUND  SPEED,  OCEAN  SURFACE  =  1500  rtfs 

SOUND  SPEED,  SOURCE  DEPTH  =  1615  rtfs 


Figure  17.  Plot  of  £(y)  versus  depth  y  for  different  launch  angles  |3o : 
n2(y)  with  negative  gradient. 
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